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ABSTRACT 


CARDINAL  LACUNARY  SPLINE  INTERPOLATION 

Following  the  ideas  of  Lipow  and  Schoenberg  on  Cardinal 
Hermite  Interpolation  (C.H. I.P. ),  we  study  here  a  class  of  Cardinal 
Lacunary  Interpolation  Problems  (C.L.I.P.).  A  special  case  for  quintic 
splines  has  recently  been  treated  by  Meir  and  Sharma  and  a  generali¬ 
zation  of  their  results  is  contained  in  the  work  of  Demko.  The  first 
Chapter  deals  in  full  with  the  statement  of  the  problem  and  a  complete 
solution  to  the  problem  of  (0 ,1, . . . ,r-2,r)  interpolation  by  g-splines. 
When  the  data  is  of  power  growth,  the  interpolation  problem  has  a 
unique  solution.  Chapter  II  is  devoted  to  the  problem  of  the  zeros  of 
the  characteristic  polynomials  which  arise  in  the  study  of  C.L.I.P. 
These  results  generalize  the  theorem  of  Lipow  and  Schoenberg  on  the 

zeros  of  II  (X)  and  supplement  it  in  some  ways.  In  Chapter  III  we 
n,r 

study  Hankel  determinants  of  Euler-Frobenius  polynomials,  their 
relation  with  the  polynomials  11^  r(A)  and  their  generalizations. 

Chapter  IV  deals  with  a  generalization  of  exponential  Euler  splines 
which  arises  in  a  natural  way  when  we  consider  C.H. I.P.  to  X  ,  X  ^  1. 

We  apply  these  results  in  Chapter  V  to  find  the  Fourier  transform  of 
B-splines  for  C.H. I.P.  introduced  by  Schoenberg  and  Sharma  in  the  case 
of  (0,1)  interpolation.  The  method  can  be  applied  to  more  general 
situations . 
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INTRODUCTION 


Recently  much  attention  has  been  given  to  interpolation  and 

approximation  by  spline  functions.  The  pioneering  monograph  of  I.  J. 

Schoenberg  [10]  of  1946,  deals  with  what  is  now  called  the  Cardinal 

interpolation  problem  (C.I.P.).  The  subject  has  been  further 

investigated  by  him  recently  in  [12]  and  [13].  Further  references  on 

C.I.P.  can  be  found  in  the  excellent  monograph  of  I.  J.  Schoenberg  [15]. 

Let  n,r  be  positive  integers  with  n  >  2r  -  1.  The  class  S  of 

—  n,r 

Cardinal  splines  of  degree  n  with  integer  knots  of  multiplicity  r 
consists  of  functions  S(x)  satisfying 

(1)  S (x)  e  Cn"r 
and 

(2)  S(x)  is  a  polynomial  of  degree  n  in  each  of  the  intervals 
[v,  v+1]  for  all  integers  v. 

In  [6]  Lipow  and  Schoenberg  and  in  [17]  Schoenberg  and  Sharma 
studied  the  Cardinal  Hermite  interpolation  problem  (C.H.I.P.). 

Problem  (C.H.I.P.).  Given  r  bi-infinite  sequences  of  numbers 

(3)  y^  =  (y^)(i  =  0,  1,  ...,  r-1)  V  integers  v, 

find  S(x)  e  S  such  that 
BjlI 

(4)  S(i)(v)  =  y^(i  =  0,  1,  ...,  r-1)  V  integers  v. 

Recently  Meir  and  Sharma  [8]  studied  the  problems  of  existence, 
uniqueness  and  convergence  of  quintic  splines  and  quintic  g-splines 
which  together  with  their  second  derivatives  interpolate  a  given 
function  and  its  second  derivative  respectively  at  a  finite  number  of 
equidistant  knots.  In  the  usual  terminology  we  shall  call  this  problem 
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the  (0,2)  interpolation  problem.  More  generally  we  shall  say  that  S(x) 

belongs  to  the  class  of  Cardinal  g-splines  of  degree  n  correspond- 

n ,  r 

ing  to  (0,  1,  r-2,  r)  interpolation  problem  when 

(5)  S (x)  e  Cn  r  1  (-00  ,00), 

(6)  S(x)  is  a  polynomial  of  degree  n  in  each  of  the  intervals 
[ v , v+1 ]  for  all  v, 

and  (7)  S(n_r+1)(v+)  =  S (n_r+1) (v-) (v  =0,  ±1,  ±2,  ...). 

In  [19]  B.  Swartz  and  R.  Varga  have  shown  that  the  error  bounds  of  Meir 
and  Sharma  are  best  possible  in  some  sense.  Using  the  technique  of 
Meir  and  Sharma,  S.  Demko  [2]  has  also  extended  their  results  to  higher 
degree  splines. 

The  main  object  of  this  investigation  is  to  study  the  problem  of 
Meir  and  Sharma  for  Cardinal  splines  of  degree  n,  using  the  approach  of 
Lipow  and  Schoenberg  [6].  However  we  do  not  go  into  the  convergence 
problem  here.  In  the  case  of  C.H.I.P.  this  approach  consists,  first  of 
all,  in  finding  the  eigensplines  of  the  space 

(8)  S  =  (S(x)  e  5  :  S^(v)  =  0  V  integers  v  (p  =  0,1, . . .  ,r-l) } . 

n ,  r  n ,  r 

If  S(x)  £  S  ,  S(x)  t  0  and  if 
n,r 

(9)  S(x+1)  =  X  S(x)  V  real  x,  X  ^  0, 

then  S(x)  is  called  an  eigenspline  (E.S.)  and  X  is  called  the  eigenvalue 

(E.V.).  These  eigensplines  are  used  in  the  construction  of  fundamental 

splines  L  (x)(s  =  0,  1,  ...,  r-1)  which  satisfy 
s 
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(10) 


(v)  =  0  V  integer  v  i  0  (p  =  0,  1,  r-1) , 

s 

< 


L(p)(0)  =  0  (p  ji  s)  and  L(s)(0)  =  1  . 

\  s  s 

The  C.H.I.P.  is  essentially  solved  once  these  fundamental  splines  have 
been  constructed. 

Here,  for  the  Cardinal  lacunary  interpolation  problems  (C.L.I.P.), 
the  processes  of  finding  the  corresponding  E.S.  and  of  constructing  the 
fundamental  splines  are  basically  the  same  as  in  C.H.I.P.  The  problem 
of  finding  the  eigensplines  reduces  to  the  study  of  the  zeros  of  some 
characteristic  polynomials.  More  precisely  let  P  =  1Kj)||  (i»j  =  0,  1, 
2,  ...)  be  the  infinite  matrix  of  the  binomial  coefficients  so  that  the 
characteristic  matrix  is 


(11) 


P  -  A  X  =  ||  Cj)  -  X  6  ||  (l,j  =  0,  1,  2,  ...). 


We  shall  denote  by  P 


/. 

iQ,  i1,  ...,  iv 

:A 

J0>  •  •  • »  3V 


the  determinant  of  the  sub¬ 


matrix  obtained  from  (11)  by  deleting  all  the  rows  and  columns  except 
those  labelled  (i^,  i^,  ...,  i^}  and  {jQ,  j^,  ...,  j^}  respectively, 


and  by  P 


v  v  \ 

j0>  3±y  •  •  •  ,  Jv 


the  corresponding  determinant  obtained  from  P. 


In  connection  with  the  C.H.I.P.  Lipow  and  Schoenberg  [6]  proved  that  the 


polynomial 


(12) 


n 

n,r 


(x) 


r+l, 

1, 


n 

n-r 


:  X 


is  a  reciprocal  polynomial  whose  zeros  are  real,  simple  and  of  sign 
(-l)r.  They  use  the  technique  of  reducing  the  problem  to  an  eigenvalue 


f? 


problem  and  applying  a  powerful  theorem  of  Gantmacher  and  Krein  on 


eigenvalues  of  oscillating  matrices.  However  this  does  not  yield  the 


interlacing  property  of  the  zeros  of  n 


(X) ,  a  phenomenon  that  occurs 


n,r 


when  r  =  1.  This  is  due  to  G.  Frobenius  (see  [14]). 

Summary  of  Thesis.  This  work  is  divided  into  five  chapters.  The 
first  chapter  gives  a  full  discussion  of  the  special  case  of 
(0,1,. . . ,r-2,r)  interpolation  by  g-splines.  The  corresponding  funda¬ 
mental  splines  are  constructed  by  a  slight  modification  of  the  method 
of  Lipow  and  Schoenberg  [6], 

In  Chapter  II  we  give  another  proof  of  the  theorem  of  Lipow  and 
Schoenberg.  Furthermore  we  show  that  the  zeros  of  II  (X)  and 


n,r 

IIn_^  r(X)  are  indeed  interlacing.  Exploiting  this  interlacing 


property  we  are  able  to  prove  the  corresponding  results  on  the 
characteristic  polynomials  arising  from  C.L.I.P. 

Chapter  III  is  devoted  to  the  study  of  the  relations  between  the 
characteristic  polynomials  of  C.L.I.P.  and  the  Euler-Frobenius  poly¬ 
nomials  II  (X)  =  II  ,  (X)  .  In  particular  the  Hankel  determinant  of 
n  n,  1 

11^ (X)  /n !  is  related  to  II ^  r(X)  of  Lipow  and  Schoenberg. 

In  Chapter  IV  we  study  the  analogue  of  the  exponential  Euler 
spline  for  the  case  of  C.H.I.P.  We  call  them  exponential  Hermite- 
Euler  splines.  The  polynomial  component  of  the  exponential  Hermite- 
Euler  spline  of  degree  n  in  the  interval  [0,1]  turns  out  to  be  a 


linear  combination  of  the  polynomials 


nential  Euler  polynomials,  see  Schoenberg  [14]).  We  investigate  the 
convergence  as  the  degree  of  the  spline  tend  to  infinity,  for  the  case 
r  =  2.  Also  we  make  a  brief  mention  of  the  analogue  of  exponential 


Euler  spline  in  relation  to  the  problem  of  (0,2)  interpolation  by  g- 
splines . 

The  exponential  Euler  splines  introduced  by  Schoenberg  [14]  are 
extremely  useful  (see  [15],  [16]).  To  demonstrate  the  usefulness  of 
the  exponential  Hermite-Euler  splines  introduced  in  Chapter  IV,  an 
application  is  made  to  compute  the  Fourier  transforms  of  the  B-splines 
and  fundamental  splines  for  C.H.I.P.  This  is  the  subject  of  Chapter 
V. 

Appendix  I  is  a  table  of  examples  of  the  characteristic  poly¬ 
nomials  in  some  special  cases.  Appendix  II  contains  examples  of  B- 
splines  for  quintic  g-splines  and  ordinary  splines  for  (0,2)  and  (0,3) 
interpolation  respectively. 


CHAPTER  I 


STATEMENT  OF  C.L.I.P.  AND  A  SPECIAL  CASE 

We  shall  begin  with  the  statement  of  the  Cardinal  Lacunary 
Interpolation  Problems  (C.L.I.P.)  to  be  considered  in  the  thesis.  In 
§2  we  state  the  main  results  for  a  special  case  whose  solution  will  be 
accomplished  in  the  course  of  §  §  3  and  4. 

1.  Statement  of  problems. 

Let  n,r  be  positive  integers  such  that  n  >_  2r-l.  Further,  let 
p,q,k  be  non-negative  integers  such  that 

(1.1)  r  =  p+q,  1  <  p  £  r,  0  <_  q  <_  r-1,  and  p  <_  k  <_  n-r-q+1  . 

We  shall  use  the  notation  k,q  =  {k,k+l, . . . ,k+q-l}  and 
0,p;k,q  =  {0,1, . . . ,p-l,k,k+l, . . . ,k+q-l} . 

s  t 

If  s,t  are  integers  with  0  <  t  <_  r  and  0  <_  s  <_  n-r+1,  let 
denote  the  class  of  Cardinal  spline  functions  S (x)  satisfying  the 
following  conditions: 

(1.2)  S(x)  is  a  polynomial  of  degree  n  in  each  of  the  intervals 
[v,v+l]  (v  =  0,±1,±2, . . .)  , 

(1.3)  S (x)  e  CS-1(-°°,°o)  , 

(1.4)  S^P^(v+)  =  S^(v-)(p  =  s+t ,s+t+l, . . .  ,n-r+t)  V  integers  v. 

If  t  =  0,  (1.3)  and  (1.4)  together  imply  that  Ss,°  =  S 

XT  ^  XT  XX  y  IT 

If  t  =  1  and  s  =  n-r  then  Sn  r’^  =  S°  ,  the  class  of  Cardinal 

n,r  n,r 

g-splines  for  (0,1, . . . ,r-2,r)  interpolation  defined  in  the  Introduction. 
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We  shall  call  the  following  problem  C.L.I.P.  of  type 


(0,p;  k,q)  in 

n,r 

PROBLEM.  Given  positive  integers  p,q,k  satisfying  (1.1)  and  r  bi¬ 
infinite  sequences  of  numbers 

0.5)  y(1)  =  (ya))  (l  e  oTp;  k7?) 

satisfying 

(1*6)  y^  =  0 ( |  v | Y)  (i  e  0  ,p ;  k,q) 

s  t 

for  some  y  >  0,  to  find  a  function  S(x)  e  S  *  such  that 

n,r 

(1.7)  S^(v)  =  (i  £  0,p;  k,q)  V  integers  v. 

Observe  that  when  k  =  p  or  q  =  0  the  C.L.I.P.  reduces  to  the 
C.H.I.P.  In  the  case  of  a  finite  interval,  instead  of  the  whole  real 
line  k,  the  corresponding  problems  of  type  (0,2)  in  S  ?  and  5^  have 
been  dealt  with  by  Meir  and  Sharma  [8],  while  the  problem  of  type 
(0,p,  k,q)  (k  =  2m-r-q)  in  -^m-l  r  ^as  ^een  treated  by  S.  Demko  [2], 

2.  A  special  case. 

We  shall  discuss  in  full  the  special  case  of  (0,1, . . . ,r-2,r) 
interpolation  by  Cardinal  g-splines.  In  our  terminology  this 

- - -  pr 

corresponds  to  C.L.I.P.  of  type  (0,r-l;  r)  in  S6  .  We  shall  prove 

n ,  r 

Theorem  2.1.  Given  r  bi-infinite  sequences  (1.5)  satisfying  (1.6)  the 
C.L.I.P.  of  type  (0,r-l;r)  has  a  unique  solution  S(x)  £  5^m_^  r  such 
that 

(2.1)  (x)  =  0 ( |  x | Y)  (p  £  0,r-l;r)  as  1  x|  ->•  00 . 


In  §4  we  shall  construct  the  fundamental  functions 


Vx) 


(2.2) 


=  L^m-l  s<x>  e  ^2m-l  r^S  e  0>h-l,r)  which  satisfy 

(  <  )  — 

L^P'(v)  =  0  (p  e  0,r-l;r),  v  =  ±1,±2,±3,...  , 
s 

.. 

L(p)(0)  =  0  (p  e  0,r-l;r),  p  /  s  and  L(s>(0)  =  1  . 
S  s 


\ 

In  terms  of  these  fundamental  functions  the  spline  S(x)  in  Theorem  2.1 
is  given  by 


(2.3)  S(x)  =  l 

v=-°° 


se  0,r-l;r 


l— 


3.  The  characteristic  polynomials. 

Following  Schoenberg,  we  set 

(3.1)  =  {S(x):S(x)  e  S ®  ,S^P^(v)=  0  (p  e  0,r-l;r)  V  integers  v}  , 

n  y  IT  n  y  IT 

0  a 

and  call  S(x)  e  S6  S(x)  t  0, an  eigenspline  (E.S.)  if  it  satisfies 

n ,  r 

the  functional  relation 

(3.2)  S(x+1)  =  AS(x)  for  all  x  e  R,  A  ^  0  . 

The  number  A  is  called  the  eigenvalue  (E.V.)  of  the  eigenspline  S(x). 

°Q 

Lemma  3.1.  The  space  ^  is  a  linear  space  of  dimension  n-2r+l. 

Proof .  It  is  clear  that  the  class  of  all  polynomials  P(x)  of  degree  n 
satisfying 

(3.3)  P(p)(l)  =  P(p)(0)  =0  (p  =  0,1, .. .  ,r-2,r) 

is  a  linear  space  of  dimension  n-2r+l.  Thus  we  have  only  to  show  that 


. 


■ 


■ 


u 

every  such  P(x)  determines  a  unique  spline  S(x)  e  Sf  To  this  end 
we  set 


n,r 


(3.4)  S(x)  =  P(x)  +  l  c^Cx-l)"  s  +  l  c''^(x-2)1'"  s;  + 


n-s 


(s) 


(n-s) 


+  I  c<f(-x)"-x  +  l  c^h-x-1)""8 


+  . 


where  the  summations  are  taken  over  all  s  z  0,r-l;r.  Then 


(p) 


(p) 


S(x)  z  S&  and  (0)  =  S  (1)  =  0  (p  e  0,r-l;r).  The  conditions 


n,r 


S^(2)  =  0  (p  e  0,r-l;r)  give  a  non-homogeneous  system  of  r  equations 

(s)  —  — - 

in  r  unknowns  c  (s  e  0,r-l;r)  whose  determinant  is  equal  to 

1!  2!  .  .  .  (r-2)  !r  !P(n  r  ,n  *  **’  *n)  which  is  non-zero.  Hence  c^ 

0  ,  1  , . . .  ,r-2 ,r  1 

(  s  ) 

are  uniquely  determined  in  terms  of  P(x).  Suppose  c  (k  =  l,2,...,v-l) 

(s') 

have  been  so  determined.  Then  c  are  uniquely  determined  in  the  same 


V 


(P) 


way  b-y  the  conditions  S  (v+1)  =  0  (p  e  0,r-l;r).  By  induction  all 
(s) 


the  c  (k  =  l,2,...,s  z  0,r-l;r)  are  uniquely  determined  in  terms  of 

K. 

(s) 


P(x)  by  the  conditions  Sv  y (v)  =  0  (v  =  2,3,...  ,s  z  0,r-l,r). 

(s) 


Similarly  the  cv  y(k  =  0 ,-1,-2 , . . . ,s  z  0,r-l;r)  are  uniquely  determined 

K. 


(p) 


by  the  conditions  Sv  (v)  =  0  (v  =  -l,-2,...,p  e  0,r-l;r). 


□ 


Eigenvalues. 

°e 

Suppose  S(x)  z  Sb  is  an  eigenspline  satisfying  (3.2).  Then 

n ,  r 

S(x)  is  uniquely  determined  by  its  polynomial  component  P(x)  in  [0,1], 
It  follows  from  (3.1)  that 


(3.5) 


,(p) 


_  P(P) 


(1)  =  P  ^'(0)  =  0  (p  e  0 ,r-l ;r)  , 


while  (3.2)  together  with  the  continuity  conditions  on  S(x)  give  rise 


to  the  additional  relations 


. 


10. 


(3.6)  P(p)(l)  =  XP(p)(0)  (p  =  r-l,r+l, . . . ,n-r ,n-r+l)  . 


By  (3.5)  we  can  write 

(3.7)  P(x)  =  a  xn  +  a  (^)x11  ^  + 

o  1  l 


n 

+  a  (  )x 
n-r-1  n-r-1 


r+1 


+  Vr+l^-r+P^'1 


Writing  equations  (3.5)  and  (3.6)  upwards  with  increasing  p,  we  obtain 
a  homogeneous  system  of  (n-r+1)  equations  in  (n-r+1)  unknowns 


a  , a, ,...,a  ,,a  , , .  The  determinant  of  the  system  is  the 

o  1  n-r-1  n-r+1  J 

(n-2r+l)-th  degree  polynomial  in  A  given  by 


(3.8)  II8  (A)  =  P(r  1,r^1’  »  n  , ,  ;  X)  = 

n,r  0  ,  1  , ... , n-r-1, n-r+1 


1 

1 

1  ( 

1 

1  ( 
1  ( 


I*'1) 

(1-X) 

0 

0  . 

(r+1) 

<“> 

(r+!) 

(1-X) . 

• 

• 

n-r-1. 

1  ; 

. (1-X) 

(n“r) 

. (  n"r 

n-r+1. 

1  ; 

✓n-r+1 
.  n-r-1 

n-r+2. 

1  ' 

, n-r+2 

(  n  ) 
Ii-r-r 


0 

0 

0 

0 

(1-A) 

n-r+2 

n-r+1 

• 

(  n  ) 

Vr+l' 


( 


) 


Thus  the  E.V.’s  of  the  E.S.'s  of  the  null  space  are  the  zeros  of 

n,r 

n8  (A) .  In  §  7  we  shall  prove  the  following 
n ,  r 


Theorem  3.2.  Let  n  =  2m-l  be  an  odd  integer  such  that  m  >  r+1 


(r  =  2,3,...).  The  polynomial  n8^_^  r(A)  is  a  reciprocal  polynomial  of 


11. 


degree  d  =  2m-2r  with  real  simple  zeros  two  of  which  are  of  sign 
(-l)r  ^  and  (d-2)  of  sign  (-l)r.  (For  proof  see  §7  p.  ). 


Let  us  label  the  zeros  of  II ^  ..  (A)  by  {A.}^  ..  so  that 

2m- 1 , r  J  J=1 


(3.9)  A,  <  -1  <  A-  <  0  <  A_  <...<  A  <  1  <  A  <  . .  .  <  A ,  . 

d  12  m-r  m-r-1  d-1 


in  case  r  is  even  (with  obvious  modifications  when  r  is  odd)  and 

( 3 . 10)  A  .  A  ,  =  A  .  A  n  =  ...  =  A  A  , ,  =  1  . 

1  d  2  d-1  m-r  m-r+1 

The  eigensplines  are  constructed  as  follows.  Let  P (x)  be  a  solution 

of  the  system  of  equations  (3.5)  and  (3.6)  when  A  =  A .  Then  set 

S . (x)  =  P . (x)  (x  e  [0,1])  and  extend  S . (x)  by  the  functional  relation 
3  3  3 


(3.11)  S.(x+1)  =  A.S.(x)  (x  e  R) . 

3  3  3 


It  follows  that  for  every  integer  v  we  have 


(3.12)  S_.  (x+v)  =  Aj Sj  (x)  (x  e  R)  . 


Lemma  3.3.  Let  S . (x)  be  an  E.S.  corresponding  to  A.. 

_ , _  .  3  -  ■ .  ,  .  .  ..I  .  ■  -  -  ■  j  ...  i  —  ■  .  -  ■  ■  ■  -■  3 

Then  S(j-1)  (0)  ^  0. 

(lT“l ) 

Proof.  Suppose  S  (0)  =  0.  Then  its  polynomial  component 


_  ,  ,  2m-l  .  ,2111-1,.  2m-2 

(3.13)  P  (x)  =  aQX  +  a^(  ^  )x  +  . 


,  2m- 1  r+1 

.  +  a  0(0  „)x 

2m-r-2  2m-r-2 


satisfies  the  relations 


(3.14)  P(^(l)  =0  (p  =  0,1,. ..,r) 


and 


■ 


12. 


(3.15)  P(^(l)  =  X^P^CO)  (p  =  r+1 , . . .  ,  2m-r-2 , 2m-r )  . 

Now  (3.14)  and  (3.15)  give  a  system  of  (2m-r)  equations  in 
(2m-r-l)  unknowns  (k  =  0 ,1 , . . . , 2m-r-2) .  Since  P^ (x)  t  0 ,  it  fol- 

(X)  and  at  the  same  time  a  zero  of 

»tTi 

r-l,r+l, . . . ,2m-2 


lows  that  X.  is  a  zero  of  H.  .  ,, 

j  2m- 1, r+1 


n  (x)  =  p 

2m-2,r 


0  ,  1  , . . . , 2m-2-r 
r 


:XJ.  But  then  (see  §7)  all  the  zeros 

r+1 


of  II 2m_2  r(^)  are  °f  sign  (-1)  except  the  zero  (-1)  ,  and  all  the 

zeros  of  II  ..  (X)  are  of  sign  (-1)  .  Since  X.  4  (-1)  these 

2m- 1, r+1  j 

(  T—  1  ) 

are  impossible.  Hence  S  \  '  (0)  ^  0. 


□ 


Remarks .  1)  The  proof  of  lemma  3.3  also  establishes  the  fact  that  the 

Or 

matrix  whose  determinant  is  II ^  ..  (X.)  has  rank  d-1.  Thus  the  eigen- 

2m-l,r  j 

splines  S . (x)  corresponding  to  the  eigenvalue  X.  are  uniquely  deter- 
1  3 

mined  up  to  a  constant  factor,  and  so  there  are  exactly  d  distinct 
eigensplines  in  r* 

2)  By  a  lemma  of  Schoenberg  and  Ziegler  ([18],  p.  424,  lemma 
2)  it  can  be  shown  that  these  eigensplines  can  have  at  most  a  finite 
number  of  zeros  in  (0,1). 

Since  the  X  (j  =  l,2,...,d)  are  distinct  the  relation  (3.12) 
gives  the  following 

Lemma  3.4.  The  eigensplines  {S.(x)}f_^  are  linearly  independent  and 

— - 5- 

so  form  a  basis  of  the  vector  space  Sr  _ 

r  2m-l,r 

We  shall  normalise  the  eigensplines  so  that 


(3.16)  S_.  (x)  >0  (x  e  (0,6))  for  some  6  >  0, 


and 


13. 


(3.17)  S(J  1) (0)  =  1  . 


Lemma  3.5.  If  S(x)  e  S? 

2m- 1 , r 


and  S(x)  =  0(|x|^)  for  some  y  >  0 


then  S(x)  =  0. 


Proof 


d 

.  By  lemma  3.4  we  can  write  S(x)  =  )  c.S.(x).  Since  I  A. 

3=1  3  3  3 


+  1, 


the  hypothesis  of  the  lemma  together  with  (3.12)  implies  that  c^  =  0 


(j  =  l,2,...,d).  Hence  S  (x)  =  0  . 


4.  Fundamental  splines  . 

Recall  that  the  fundamental  splines  L  (x)  E 

s 

L?  —  (x)  e  (s  e  (0,r-l;r))  should  satisfy  the  following 

2m- 1 , s  n ,  r 


conditions 

/ 


(4.1) 


L(p) 

s 

u(p) 


(v) 

(0) 


=  0 

=  0 


(p  e  0 ,r-l;r) , 


(p  e  0,r-l;r , 


\ 


v  =  ±1,±2,  . . .  , 

p  ^  s)  and  L (0) 
s 


1  . 


We  shall  construct  these  fundamental  functions  in  terms  of  the  eigen- 
splines  obtained  in  §3.  For  this  purpose  let  us  consider  the  following 
functions 


(4.2) 

F1(x) 

m-r 

=  y  c.s.(x)  , 

j=l  3  3 

(4.3) 

F2(x) 

2m- 2  r 

=  l  d  S  (x) 

j=m-r+l  J  J 

(4.4) 

P(x) 

2m- 1 

2m- 2 

=  a  x 
o 

+  a^x 

) 


+  . . .+ 


a0  0x 
2m-r-2 


r+1 


+  a 


r-1 

2m-rX 


(s  e  0,r-l;r)  , 


and 


14. 


//  c\  n /  \  ,  2m-l  2ra-2  r+1  r-1  ,  x 

(4.5)  Q(x)  =  b  x  +  b  x  +  . .  .+  b~  _x  +  b„  x  +  — 

o  1  2m-r-2  2m-r  s 


(s  e  0,r-l;r) 


Observe  that  F.(x)  e  Sf  -  (i  = 
i  2  m- 1  y  it 

(x)  e  L  (-«>,-l].  Also  P^(0) 
and  P(s)(0)  =  Q(s)(0)  =  1. 


1,2)  and  F^(x)  e  L  [1,°°)  while 
(p) 


=  Q  (0)  =  0  (p  e  0,r-l;r  and  p  ^  s) 


If  we  set 


(4.6) 

P(p)(0)  = 

q(p> 

(4.7) 

•O 

/— \ 

H* 

'w' 

II 

m-r 

I 

and 

Q(P)  (-1)  ■ 

j-i 

2m- 

(4.8) 

-  I 

:.S^ 
3  3 


j=m-r+l 


d.S(p)(-l)  (p  =  0,1,... ,2m-r-2,2m-r) 
3  3 


we  obtain  a  non-homogeneous  system  of  6m-4r  equations  in  6m-4r  unknowns 

{ a . ,b . : j =0 , 1 , . . . , 2m-r-2 , 2m-r}  ,  {  c . } m  *  and  {  d . } .  We  want  to 
33  3  1=1  3  j=m-r+l 

show  that  this  system  is  non-singular.  Observe  that  the  constant  terms 

1  (-1)S_P 

of  this  system  of  equations  consist  of  — - ry  ,  - r-y —  which  appear  in 

the  last  two  sets  of  equations  (4.7)  and  (4.8)  for  p  =  0,1,. ..,s,  and 

the  rest  are  zeros.  The  corresponding  homogeneous  system  corresponds  to 

the  system  of  equations  given  by  (4.6),  (4.7)  and  (4.8)  in  which  the 
s 

terms  — j-  are  absent  in  P(x)  and  Q(x).  If  this  system  is  singular  we 

S  • 

°e 

would  obtain  a  spline  S(x)  e  SB  ,  S(x)  t  0,  and  S(x)  e  L  (-00,00) .  This 

XI  y  XT  _L 

is  impossible  because  of  lemma  3.5.  Hence  the  above  system  is  non- 

a 

singular.  This  means  that  we  can  find  a  unique  spline  L  (x)  e 

s  zm— ± , r 

(s  e  0,r-l;r)  satisfying  (4.1).  Furthermore  L  (x)  e  L  . 

S  JL 


The  symmetry  relation 


(4.9)  L  (x)  =  (-1)* * * * S  L  (-x)  (s  e  0,r-l;r) 

s  s 

follows  from  the  following  lemma  whose  proof  is  immediate. 

Lemma  4.1.  Let  L  (x)  =  L  (-x) .  Then  L  (x)  e  and  satisfies 

 s  s_ s_ n,r 

L^Cv)  =  0  (p  e  0,r-l;r)  V  integers  v  ^  0,  L^P^  (0)  =  0  (p  e  0,r-l;r, 
s  s 

p  V  S)  and  L(s)(0)  =  (-1)S  . 

s 

Lemma  4.2.  There  exists  positive  constants  A  and  a  depending  on  m  and 
r  such  that 

(4.10)  |L^p\x)  I  <  A  e  (p  e0,r-l;r)  . 

Proof.  Recall  that  L  (x)  is  constructed  so  that 
-  s 

m-r 

(4.11)  L  (x)  =  ][  c.S.(x)  (x  >  1)  . 

S  L  i  i  — 

j=l 


Hence,  in  view  of  the  relation  (4.9),  it  is  clearly  enough  to  prove  the 

same  estimate  for  S^P^(x)  (j  =  1,2,..., m-r)  for  x  >_  0. 

3 

Suppose  v  <_  x  <  v+1.  By  (3.12)  we  have 

(4.12)  S(p)(x)  =  XV.  S(p)(x-v) 

3  3  3 


(P) 


v 


and  so  S  .  (x)  <  BA.  <  B|A. 


-li 


(j  =  1,2 , . . . ,m-r) ,  where 


B  =  sup  S^P^(x).  Hence 
0<x<l  J 


(4.13)  | S(p) (x) |  £  B 


ex 


p{  -  x  log | A . j  h 


A 


(j  =  1,2,. . . ,m-r) . 


If  we  set  a  = 


mm 


log | A . |  \  we  have 


j-1,2, . . . ,m  r 


16'. 


(4.14)  |  (x)  |  <_  B  — j1— •  (j  =  1,2 , . . .  ,m-r)  ,  x  >_  0  . 


Hence  (4.10)  follows  from  (4.11)  and  (4.14). 


□ 


Proof  of  Theorem  2. 1. 


We  shall  show  that  if 


(P) 


(4.15)  =  0  ( |  v  | Y)  (p  e  0,r-l;r)  for  some  y  >  0, 


v 


then 


(4.16)  S(p)(x)  -l  l 


s  e  0,r-l,r 


y(s>L(p)(x-v) 

J  V  S 


(p  e  0 ,r-l;r) 


converges  locally  uniformly  and 


(P) 


(4.17)  S VH/  (x)  =  0(|x|Y)  (p  e  0  ,r-l;r)  as  |x|  -> 


The  proof  of  (4.16)  and  (4.17)  is  the  same  as  in  [13]  p.  415 
with  obvious  modifications.  □ 


CHAPTER  II 


PROOF  OF  THEOREM  3.2  AND  SOME  GENERALISATIONS 

p 

In  §3  the  zeros  of  E*  (X)  are  used  in  the  construction  of 

2m-l,r 

the  eigensplines  for  the  C.L.I.P.  of  type  (0,r-l:r)  in  ,  .  This 

2m-l,r 

chapter  is  devoted  to  the  proof  of  theorem  3.2.  We  shall  also  discuss 
the  corresponding  results  for  some  general  cases.  These  results 
generalise  the  theorem  of  Lipow  and  Schoenberg  [6]  and  also  our  theorem 

3.2. 


5.  The  characteristic  polynomials  for  general  C.L.I.P.. 

Let  n  be  a  positive  integer  such  that  n  >_  2r-l  and  p,q,k  be 
non-negative  integers  satisfying  (1.1).  Further,  let  s,t  be  non¬ 
negative  integers  with 

(5.1)  k+q  s  <_  n-r+1  ,  0  t  <  r  . 

Following  the  same  argument  as  in  §3,  we  see  that  for  the  C.L.I.P.  of 

type  (0,p;  k,q)  in  S  *  the  E.V.'s  of  the  E.S.’s  belonging  to  the  null 

n ,  r 

space 

(5.2)  Ss’l  (07p;  k^q)  =  ( S (x)  e  S®**  :  S(p)(v)  =  0 


for  p  e  0,p;  k,q  and  v  =  0,±1,±2,...} 


s  t  ~  ■  — — 

are  the  zeros  of  the  characteristic  polynomials  n  ’  (0,p;  k,q:  X)  where 

n ,  r 

„  s ,  t „  ,r-t ,r-t+l , . . . ,n-s-t ,n-s+l ,  ...  ,  n  .  N 

(5.3)  nn;r(0,p;  k,q:x)  =  P(  0  ;  x  ;n_k_q;n_k+1 . 


17 


_ 


18. 


n-s-r+l< 


s-k-q 


q< 


k-p^ 


(r^> . . .  (l-X)  0 


.0 


1  (n  ®  t) .  (1-X)  0 


1  (n-k-q+1) 


1  rf1) 


(?) 


,0 


1  (n  l+1) . (1-X)  0 


.0 


(n  k  q)  .  (1-X) 


0 


,0 


n-k-q+1  Q 
n-k-q 


0 


(n  k+i)  (1_A)  Q  _Q 

n-k-q 


rn  P')  .  .  (  n  p  'i  r  n  p  ^  n-n 

^  1  ;  .  n-k-q'  Vk+r,,U  ; 


n 


C  n  )  ,  n  )  (  n  )  (  n  )  (  n  )  (  n  ) 

*  *  V-t;  *  *  n-s-r  *  ’  n-s+1'*  ’''n-k-q'  ^n-k+l'  *  *  ; 

_ _ a - v - '  s - 


£zp^ 


r-t 


n-s-r+1 


s-k-q 


k-  p 


We  shall  write 


r 


(5.3a)  + 


rs  ,o 


n  *  (0,p;  k,q:X)  =  n  (0 , p ;  k,q:X) 


n,r 


n,r 


ns,t(0,p;  k,0;X)  =  nf^(X)  . 
n ,  r 


n,r 


Observe  that  X  occurs  in  (n-2r+l) rows  and  (n-2r+l) columns , 


rs>t 


and  it  is  clear  that  II  ’  (0 , p ;  k,q:X)  is  a  polynomial  in  X  of  degree 

n ,  r 

d  =  n-2r+l.  By  the  same  reason  as  for  n  (X) ,  the  polynomials 


n,r 


g  |-  ■  —  —  ■  -  — ■  - 

II  *  (0,p;  k,q:X)  are  reciprocal  polynomials.  Some  particular  cases  are 
n ,  r 


. .  -  T 


•* 
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(  TTS>° 


n  ’  (0,p;  k,0:X)  =  n  (X)  , 


n,r 


n,r 


(5.4)  < 


Jln  r,1(0,r-l,r:A)  =  II8  (A)  . 


\  n,r 


n,r 


For  the  sake  of  convenience  we  shall  also  write 


■p  __  y  1  _  ^  ^ 

(5.6)  n_  _»  (0,p;k,0:A)  =  n  (A) 

n,r 


n,r 


and 


(5.7)  IIS,O(0,r-l;r:A)  =  II  (A)  . 
"  ”  n,r 


n,r 


The  relation 


*S»  ^ 

(5.8)  (n-r+l)n  (A)  =  r  n  (A) 


n,r 


n,r 


is  a  particular  case  of  the  following  lemma, 


Lemma  5.1.  Let  n  >.  2r-l  and  p^kjS^t,  satisfy  (1.1)  and  (5.1).  Then 


(5.9)  ns,t(0,p;k,q:A)  =  CE 


n,r 


n-q-k+1 , q 
n,r 


(0 ,r-t ,n-s-t+l, t : A)  , 


where  C 


( n )(  n .)...(  ?+1)<n.  n)“-00 

n-p  n-p-1 _ n-k+1  n-k-q  1  0 

(  n  )  (  n  )  (  n  )  (  n  )  (n)  fn) 

ln-r+tj  n-r+t-r  *  *  *  ''s+t'  ^s-l}  ’  *  *  V  0' 


Proof.  By  taking  transpose  and  the  interchanging  rows  and  columns  in 


s  t  -  — - 

the  determinant  (5.3)  we  can  write  n  ’  (0,p;k,q:A)  = 

n  >  r 


20. 


(  n  )  (n  1) 

n-p  n-p 


<  n  >  <  n_1  > 

vn-p-l;  vn-p-l' 


<  n  )  <  n'1  ) 
Kn-k+l]  ln-k+l; 


(  n  )  (  n  ^  ) 
^n-k-q'  ^n-k-q' 


(?)  O 


1 


(1-X)  0 


(1-X) 


0 


0 


0 


o 


( 


0 

(1-X) 

r-t+1 

r-t-1 

r-t+1 


0 


0 


0 

0 

)  (1-X) 

r-t 


^r-t-1^  ^ r-t-1 ^ 


,n-s+L  ,n-s-t.  ,r-t+l.  ,r-t. 

1  '  1  /•••V  -j^ 


If  we  factor  out  (JM  ,  Cp-1} . ^-k+d  >  Ck-q5 . (P  »  0  £r0m 

the  1st,  2nd  and  so  forth  until  the  last  row  respectively,  and  apply  the 

relation 

/k\  /  /Hi\  /Til” t  y  *  m  .  , 

V  /(d  -  W/  W  <”ik> - 

then  factor  out  1/ (n_"+t)  (n_"+t-l>  •  •  •  <-s+l+1>  (s+t)  Cd  • '  •  <1>  0  fr°m  the 


columns  we  obtain  (5.9). 


□ 


6.  A  determinantal  identity. 

For  convenience  of  reference  we  now  state  a  determinantal 
identity  and  a  theorem  of  S.  Karlin  which  will  be  used  in  the  sequel. 


- 


4 


Consider  the  following  n  x  (n+2)  matrix 


(6.1) 


ai  bi 


a2  b2 


(1) 

1 

f(2) 

1 

f(n-2) 

•  •  •  I  2 

(1) 

f<2> 

f (n-2) 

•  •  •  I  ^ 

abed  f(1)  f(2)  ...  f^n  2) 

n  n  n  n  n  n  n 


Let  a,  b,  c,  d,  f  f^n  ^  denote  the  column  vectors  of  (6.1) 


and  D(a,b,f)  =  D(a,b,f 


(1) 


(n_2) 

fv  )  denote  the  minors  of  (6.1) 


formed  from  the  columns  a,b,f^,  . ..,  f^n  2\  Then  an  identity  of 
S.  Karlin  ([4],  p.  7)  states  that 

D(a,c,f)  D(b,c,f) 


(6.2) 


D(a,d,f)  D(b,d,f) 


=  D(a,b ,f )D(c,d,f) 


Theorem  6.1  ([4]  ,  p.  85) .  Let  A  =  ||  a_  (|  (i,j  =  0,1, .  . .  ,n)  be  a  triangu¬ 
lar  matrix,  i.e.  a. .=0  for  j  >  i.  If  A(^ * *  *  *  * ’ i+P)  >  0  for 

ij  0 ,  1  , . . .  ,  p 


0  £  i  £  n-p,  0  <_  p  <_  n,  then  A  is  totally  positive.  Furthermore 


22. 


VV*  •  *  * ±v 
•  •  »3V 


>  0  for  0  < 


i  <  i.  <  . . .  <  i 
o  1  v 

j  <  j  j 

J  O  J1  J  V 


<  n 


whenever  j  <_  1  (k  =  0,1, ...  ,v;  v  =  0,1,. . . ,n) . 

In  the  case  of  the  matrix  P  =  II  (j)  !l  (i,j  =  0,1,2,. 


.),  it  is 


easy  to  show  by  induction  that  (see  [6],  p.  284) 

(i,p  =  0,1,2,. ..) . 


(6.3)  p(i,i+l,...,i+P)  =  x 

(J,  i  ,...,  p 


Consequently  theorem  6.1  gives  the  following 

Lemma  6.2.  The  matrix  P  =  ||(^)||  (i,j  =  0,1,2,...)  is  totally  positive. 
Furthermore 


(6.4)  P 


W**‘,:Lv 
Vjl’**‘  ,jv 


>  0  for  0  < 


i  <  i.  <  . . .  <  i 
o  1  v 

i  <  j  _<...<  j 
Jo  J1  Jv 


whenever  j  <_  i^  (k  =  0,1, . .  .  ,v)  . 


We  deduce  from  (5.3)  and  lemma  6.2  that 


(6.5)  nS,t(0,p;  k,q;0)  >  0  . 


n,r 


7.  Proof  of  Theorem  3.2  . 

First  we  establish  some  identities  for  the  characteristic 
polynomials . 

Lemma  7.1.  Let  r  =  1,2,...  .  For  n  >_  2r+l  we  have 

(7.1)  m  J_1(x)n  0  Ax)  =  n{n  .  (x)}2  -  (n-r)n  .  (x)n  (A) 

n,r+l  n-2,r-l  n-l,r  n-2,r  n,r 


Let  r  =  2,3,...  .  For  n  ^  2r+l 

(7.2)  n  1(x)n  _  (X)  =  n*  (x)n  . (x)  -  n  (x)n*  .  .  (x) 

n,r-l  n-l,r  n,r  n-l,r-l  n,r  n-l,r-l 


23. 


and 


^  -y 

(7'3)  "n  r-l(X)nn  r+l(X)  =  L  ,-<X)It  rM  ~  11  ,<X)n  JX~> 

^  j  ^  -L  n  >  IT  1  J-  II  j  IT  H  $  ^  H  j  IT  n.  j  IT 


Proof .  To  prove  (7.1)  let  us  consider  the  following  (n-r+1)  x  (n-r+1) 
matrix. 


(7.4) 


1 

(*)  ... 

,  r 
(r-l) 

(1-A) 

0 

. . .  0 

0 

I 

i 

• 

• 

(rtX)  ... 

r+i 

vr-r 

(r+i) 

(1-A) 

• 

. . .  0 

• 

• 

• 

0 

. 

l 

ri’h . . . 

‘“S1) 

(n-r-!) 

n-T-1 
k  r+1  ; 

• 

• 

• 

...  (1-A) 

0 

l 

• 

(n~r)  ... 

CP 

• 

(n-r) 

• 

n-r 

vr+l; 

• 

n-r 

*  *  *  Vr-1; 

• 

(1-A) 

• 

• 

• 

1 

i 

(*)  ... 

• 

• 

(n) 

r 

• 

• 

(  n  ) 
r+l' 

• 

• 

...(  n  J 

n-r-1 

o 

Let  us  denote  the  first  and  the  last  columns  of  (7.4)  respec¬ 
tively  by  c  and  d  and  its  v-th  column  by  f^V  ^  for  v  =  2, 3,..., n-r. 

T  T 

Set  a  =  (1,0,..., 0)  and  b  =  (0,...,0,1)  .  Then  after  some  simplifica¬ 
tions  we  have  D(a,b,f)  =  (-l)n  V  ^(n  ^)II  0  (A),  D(a,c,f)  =  II  (A) , 

\  „  J  r  n-2,r  n,r+l  * 

D(a,d,f)  =  (-l)n"r“1(^)n  (A),  D (b , c , f )  =  (-l)n~ril  (A), 

D (b  , d ,  f )  =  -(““hn  _9  ,  (A)  and  D(c,d,f)  =  (-l)n'r_1n  (A)  . 

w  n,  L“  _L  n  Z  jL  -L  <w  a>  XI  y  XT 

Substituting  these  into  (6.2)  we  obtain  (7.1). 

T 

In  order  to  prove  (7.2)  we  again  set  a  =  (1,0,..., 0)  , 

T  (v) 

b  =  (0 ,0 , . . . ,0 ,1)  and  let  f  denote  the  v-th  column  of  (7.4)  for 


v  =  1,2, . . . ,n-r-l,  with  c,  d  representing  the  last  two  columns.  Then 


24. 


D(a,b,f)  = 

D(a,d,f)  = 
* 


_IIn  1 

n-1 ,  r 


“  (-1)n~r~lnn,r+l(X)’ 

(-l)n"r_1n*>r+1(X),  D(b,c,f)  =  -nn.1>r(^)»  D(b,d,f)  = 

and  D(c,d,f)  =  II  (X)  . 

~  ~  ~  n,r 


Replacing  r  by  (r-1)  in  the  above  determinants,  and  using  the 
identity  (6.2)  we  obtain  (7.2). 

In  order  to  prove  (7.3)  we  use  the  same  vectors  as  in  the 

T 

proof  of  (7.2)  except  that  b  is  replaced  by  the  vector  (0 ,1,0, . . . ,0)  . 
Then  (7.3)  is  obtained  by  replacing  r  by  r-1.  □ 


Lemma  7.2.  If  n  ^  2r-l,  then 

(7.5)  n  (A)  >  0  for  (-l)r+1A  >  0  . 
n,r  — 


Proof .  A  straightforward  computation  shows  that 


(7.6)  n_  _(A)  =  a  [(-l)r+1A]n  2r+1  +  a  [(-l)r+1A]n  2r  +. . .+  a  ,  , 

o  1  n-2r+l 


n,r 


where 

(7.7) 


r+v-jr+v-, — , r+v  ,n-r+l , . . . ,  n 

^  __  y  p  /  i  ^  ^  \ 

k  ^  0,1  , — ,r-l,r+v^,r+v2, . . . ,r+v 


(k  =  1, 2, . .  .  ,n-2r+l)  the  summation  is  over  all  the  (n  choices  of 


(V1,V2  *  *  *  * ,Vk)  ^rom  • • • ,n-2r} ,  and 


(7.8) 


a  =  p/n-r+1,.... 
o  k  0  ,  1,. . .  ,r-l 


’  11  ) 


Hence  (7.5)  follows  from  lemma  6.2. 


□ 


Theorem  7.3.  If  r  >  1,  then  for  n  >  2r+l  the  zeros  II  (A)  are  real, 
 ’ n ,  r  


simple,  of  sign  (-1)  and  interlace  with  the  zeros  of  II  _  (A). 

n-i  ^  it 


25. 


Proof .  The  proof  follows  by  induction  on  n.  We  shall  give  the  proof 

only  for  the  case  when  r  is  even.  The  case  when  r  is  odd  can  be 

treated  in  the  same  way.  We  shall  denote  the  zeros  of  II  (X)  by 

n,r  J  1 

If  we  set  n  =  2r+l  in  (7.1)  and  take  into  account  that 
n2r+l,r+l(X)  E  n2r-l,r(X)  E  1>  we  obtain 

(7.9)  (r+Dn2r+1>r(A)  -  (2r+l){n2r;r(A)}2  -  m  ^(A)  . 


N°wn2r  (A)  is  a  reciprocal  polynomial  of  degree  1  which  is  positive  for 

X  <  0.  Hence  II  (1)  =  0,  and  it  follows  that  (r-hl ) II (1)  = 

—  2r,r  '  2r+l,r 

-rll  (1)  <  0  by  lemma  7.2.  But  then  II  (0)  >  0  and  therefore 

zr-±,r-i  Zr+l,r 

^2r+l  ^aS  a  zero  (0,1).  Since  it  is  a  reciprocal  polynomial  it 

also  has  a  zero  in  (1,°°).  Thus  the  zeros  of  II  (X)  interlace  with 

z  it “r  j.  ,  r 


the  zero  of  II „  (X) . 

2r  ,r 


Let  us  suppose  that  the  assertion  has  been  proved  for  n^_2 

and  nn_x  r(X).  More  precisely,  let  us  suppose  that  the  zeros  of 

II  _  (X)  and  II  ..  (X)  satisfy  the  relation 

n-2,r  n-l,r 

(7.10)  0  <  X(a_1)<  X(a_2)<  X(a_1)<...<  X(a_2)<  X(a"^  (d  =  n-2r+l). 

112  d-Z  d-1 

It  follows  from  (7.1)  that 

(7.11)  -(n-r)n  (A  (n~1))n  (A  <n_1))  = 

IT  y  L.  -L  II  fc-  J  i  jL 

rn  n.r+l(Al(n'1>)\-2,r-l(Al(n'1))  (1  =  1«2 . ^ 


(A) 


Since  r  is  even,  lemma  7.2  says  that  the  R.H.S.  of  (7.11)  is  positive  for 


i  =  1,2, .. . , d-1.  Hence 

(n-1) 


(7.12)  n  (X. 

n,r  x 


)n  _  (X.(n  1})<  0  (i  =  1,2, .. . ,d-l)  . 

n-z  ,r  l 


26. 


It  follows  from  (7.10)  that  II  (X.^n  have  alternating  signs  for 

n-r  1 


i  =  l,2,...,d-l.  Hence  II  (X)  has  at  least  one  zero  in  each  of  the 

n ,  r 

intervals  (X^  \x^^^)(i  =  1,2, . . .  ,d-2) .  By  inductive  hypothesis 

and  the  fact  that  II  0  (0)  >  0,  we  have  II  .  (X.  >0.  It 

n-Z  ,r  n-Z ,r  1 

follows  that  II  (X  <n  ^)<  0.  But  then  II  (0)  >0.  Hence  II  (X)  has 
n,r  1  n,r  n,r 

a  zero  in  (0,X^n  ^)  and  since  it  is  a  reciprocal  it  must  also  have  a 


(n_i) 

zero  in  (Xv^  ^  ,  00 ) .  The  assertion  is  proved  by  induction. 


□ 


* 


Theorem  7.4.  If  r  >_  2,  then  for  n  >_  2r+l,  II  (X ) (=t II  (0,r-l;r:X))  has 


njLr 


r-l 


real  simple  zeros,  one  of  which  is  (-1)  The  remaining  d-1  zeros 


are  of  sign  (-1)  and  interlace  with  the  zeros  of  II  (X)  . 


ni r 


Proof.  We  shall  again  assume  that  r  is  even,  and  let  us  suppose  that 

the  zeros  of  II  ..  (X)  and  II  (X)  are  given  by 

n-i ,r  n,r 

(7.13)  0  <  X(^  <  X(*_1)<...<  X(”:J}  <  X(*}. 

We  can  certainly  assume  this  by  theorem  7.3.  Now  (7.2)  of  lemma  7.1 


implies  that 
(7.14)  n 


(x(n))n  (x(n)) 

n,r-l  l  n-l,r  l 


n*  a(n))n 


n,r  '  i  n-1 , r-l '  i 


(i  =  1,2, . . . ,d) . 


Since  II  ,  (X^)  and  II  1  ,  (X^)  are  positive  for  all  i  =  l,...,d, 

n,r-l  i  n-1, r-l  l 

it  follows  from  (7.14)  that 

(7.15)  Sgn  II  )  =  Sgn  . 


X 

By  the  same  argument  as  in  the  proof  of  theorem  7.3,  II  (X) 

n ,  r 


' 


00 


) .  For 


has  exactly  one  zero  in  each  of  the  intervals  (A^n\  if^), 

(i  =  1,2, . . .  ,d-l)  .  It  cannot  have  a  zero  in  (0,1^)  or  (A^\ 

if  it  has  a  zero  in  (0,  A^^)  it  must  also  have  a  zero  in  (A^^,00  )  and 

1  d 

vice  versa,  which  is  impossible  since  it  is  only  of  degree  d.  Since 


II  (A)  is  a  reciprocal  equation  the  remaining  zero  must  be  -1.  The 
n ,  r 

theorem  is  thus  established. 


□ 


Theorem  7.5.  If  r  >_  2,  then  for  n  >_  2r+l,  II8  ^(A)  (=  IIn  r ’ ^(0,r-l;r ; A)) 

rl  r"  ,  T 


njL.r 


r — l  x 

has  real  zeros  two  of  which  are  of  sign  (-1)  and  (d-2)  of  sign  (-1)  . 


If  n  is  odd  the  zeros  are  simple  and  are  separated  by  the  zeros  of 

"k  —  1 

IT  (A)  .  If  n  is  even  (-1)  is  a  double  zero  and  the  remaining  (d-2) 


n,r 


zeros  are  simple  and  interlace  with  the  (d-1)  zeros  of  II  (A)  of  sign 


ni.r 


(-1) 


r . 


Proof.  Again,  we  shall  assume  that  r  is  even. 


(i)  n  is  odd 


/  \  ^  (r0 

Let  us  denote  by  (p  with  P  ^  =  the  zeros  of 


II  (A).  By  theorem  7.4,  we  can  write 
n,r 

(7.16)  u(f  <l<f  <  y(^<...<  „<“>  <  X(“>. 


From  (7.3)  of  lemma  7.1,  we  have 


(7.17)  n8  (p(n))n  (p(n))  =  -n_  „  .(p^n 


n,r  "  i  n,r  l 


n,r-l  i  n,r+l  i 


(p(n))(i=l,2,...,d). 


Since  p^  >  0  for  i  =  2,3,...,d,  it  follows,  using  the  usual  argument, 
i 

that  II8  _ (A)  has  at  least  one  zero  in  each  of  the  intervals  (p^,pf^). 


n,r 


(i  =  2,3, . . .  ,d-l)  .  Thus  II8  (A)  has  at  least  (d-2)  distinct  positive 

n ,  r 


zeros . 


■ 


T 


cr 

We  shall  show  that  II ^  ^(X)  possesses  two  distinct  negative 
zeros.  Let  us  observe,  first  of  all,  that  since  II  ,(0)  and 

n  y  1.  I  JL 


-1)  = 


II  _  (0)  are  both  positive,  it  follows  that  Sgn  II  ( 

n ,  r  x  n ,  rd*  1 

(_D  (n-2r-D/2  ^  Sgn  n  _  (_i)  =  (-D  (n-2r+3)/2^  Hence>  from  (7>17) 

n  y  IT  ““  _L 

and  the  fact  that  II  (-1)  is  positive,  it  follows  that  II®  (-1)  is 

n ,r  n  ,r 

negative.  But  then  IT®  r(0)  is  positive,  and  therefore  II®  ^(X)  must 
have  a  zero  in  (-1,0).  Since  it  is  a  reciprocal  polynomial  it  must 

<y 

also  have  a  zero  in  (-°°,-l).  Thus  we  have  shown  that  II6  (X)  has  at 

n ,  r 

least  two  distinct  negative  zeros  and  at  least  ( d— 2)  distinct  positive 
zeros.  Since  its  degree  is  d  we  conclude  that  it  has  exactly  two 
distinct  negative  and  (d— 2)  distinct  positive  zeros  which  interlace 


with  the  zeros  of  II  (X)  . 


n,r 


(ii)  n  is  even 

In  this  case  II  n  (-1)  =  II  ,  .  (-1)  =  II  (-1)  =  0.  Since 
n,r-l  n,r+l  n,r 

II  (-1)  is  positive,  it  follows  from  (7.3)  that  (-1)  is  a  zero  of 
n ,  r 

IIg  (X)  .  By  the  same  argument  as  case  (i)  we  see  that  II®  (X)  has 
n,r  J  n,r 

exactly  (d-2)  distinct  positive  zeros  which  interlace  with  the  (d-1) 

•k  Q 

positive  zeros  of  II  (X).  Since  II  (X)  is  a  reciprocal  polynomial 
^  n,r  n,r 

the  remaining  zero  must  be  -1.  Thus  -1  is  a  zero  of  multiplicity  2. 
Hence  theorem  7.5  is  established.  □ 

Theorem  3.2  follows  easily  from  theorem  7.5. 


8.  Some  generalisations 

In  this  section  we  shall  discuss  some  generalisations  of  the 
results  in  §7.  More  precisely  we  shall  discuss  the  zeros  of  the 


. 


29. 


s  t  -  - 

reciprocal  polynomials  n  9  (0,p;  k,q:x)  under  certain  restrictions  on 

n ,  r 

the  integers  k,  q  and  t.  These  results  can  be  accomplished  by  using 
the  same  techniques  employed  in  §7. 

The  following  theorems  generalise  theorem  7.3. 


Theorem  8.1  If  r  >  1,  then  for  n  >  2r+l  the  zeros  of  II  (O.p: 

 —  — n,r 


n-r-q+l,q:A)  are  real, 


zeros  of  n  +1<0,p; 


simple,  of  sign  (-1)^  and  interlace  with  the 


n-r-q,q+l:A)  (respectively  II  (0,p;n-r-q,q:A)). 

n—  -L ,  r 


Theorem  8.2  Let  p,  q,  k  be  integers  satisfying  (1.1),  s,t  be  integers 

with  0  <  t  <  r  <  k-fq  <  s  <  n-r  and  suppose  that  q,t  are  both  even. 

Then  the  zeros  of  the  polynomials  II  ’  (0,p;  k,q:A),  II  ’  (0,p;  k,q:X) 

n «  it  XT"~  1  *  r 


and  nS_|5t(0,p;  k-l,q:A)  are  real,  simple  and  of  sign  (-1)^.  Further- 
n~ L  j  XT 

s  t  -  - 

more  the  zeros  of  H  ’  (0,p;  k,q:A)  interlace  with  the  zeros  of 


rs,t 


s-i,  t 


II“’“  (0,p;  k,q:A)  (respectively  II ^’"(Ojp;  k-l,q:X)). 

n-1 ,r  n-l,r 


When  q  =  0,  theorem  8.1  not  only  reduces  to  theorem  7.3  but 
also  says  that  the  zeros  of  II ^  ^(X)  and  II  r+^(0,r;  n-r,l;X)  are 
interlacing.  Similarly  when  t  =  q  =  0,  theorem  8.2  gives  theorem  7.3. 

Theorems  8.1  and  8.2  are  proved  in  the  same  way  as  theorem  7.3 
using  respectively  the  identities 


(8.1)  (n+l-p)IIn+1  r+1(0,p;n-r-q,q+l:X)nn_1^r(0,p;n-r-q,q;X) 

-  (n+l)nn  r+1(0,p;n-r-q-l,q+l:X)nn  r(0 ,p;n-r-q+l,q:X) 

=  (r+q+l)nn_1  r  (0,p-l;n-r-q-l,q+l:X)  r+1(0  ,p+l;n-r-q+l,q:X) 

(n  >  2r+l,  r  =  1,2,...)  , 


' 


30. 


and 

(8.2) 


(n-p)n®^(0,p;k,q;X)nf  ^(0,p;k-l,q:A) 


n,r 


=  nn 


s,t 


n-2  ,r 


(0  ,p;k,q:  A)IIS  t(0,p;k-l,q: 
n-l,r  n-l,r 


X) 


-  (r-t)n^^+1(O>p+l;k!qa)n®_^_1(0,p-l;k-l,q:X) 


(0  <_  t  <  r  <  k+q  s  <_  n-r)  . 


Identities  (8.1)  and  (8.2)  can  be  derived  from  (6.2) 


Theorem  7.4  admits  the  following  generalisation  which  can  be 


proved  in  the  same  way. 


Theorem  8.3  Let  n,  r,  s,  t  be  integers  with  0  <  t  <  r-1  <  s-2  <  n-r-2 , 

rs,t 


and  suppose  t  is  even.  Then  II  *  (0,r-l;r:A)  has  real  simple  zeros,  one 


nar 


r-1  f 

of  which  is  (-1)  The  remaining  (d-1)  zeros  are  of  sign  (-1)  and 


s  t 

interlace  with  the  zeros  of  II  ’  (A)  . 

n,r 


Using  theorem  8.3  we  shall  prove  a  more  general  theorem,  viz. 


Theorem  8.4  Let  n,  r,  k,  s,  t  be  integers  with 
0  <  t  < r-1  <  k  <  s— 1  <  n-r-1,  and  suppose  that  t  is  even.  Then 
II  ’  (0,r-l;k:A)  has  real  simple  zeros,  (k-r+1)  of  which  are  of  sign 


n,r 


1T“  IT 

(-1) _ and  (d-k+r-1)  of  sign  (-1)  .  Furthermore  the  zeros  of 

S  ,t:  (0,r-l;k:  A)  and  IIS  * t  (0  ,r-l;k+l :  A)  are  separated  for  k=r-l,r , . . .  ,s-2. 


n 


n,r 


n,r 


Before  proving  the  theorem,  we  first  establish  the  following 

Lemma  8 . 5  Let  n,  r,  k,  s,  t  be  integers  such  that 

0  <  t  <  r-1  <  k-1  <  s-3  <n-r-3.  Then 


- 


31. 


(8.3)  nS,t(0,r-l;k:A)nS  5-,t:(0,r-2;k-2;k:A) 
n,r  n-l,r 


-  nS,t(0,r-l;k+l:A)nS“^,t:(0,r-2;k-2,2:A) 
n,r  n-l,r 


=  nS,t(0,r-l;k-l:A)nS  (0 ,r-2 ;k-l  ,2  :  A) 

n,r  n-l,r  ’ 


(8.4)  IIS,t(0,r-l;k:A)nS,t  (0  ,r-l;k-l;k+l:  A) 
n,r  n,r+l 


-  nS,t(0,r-l;k+l:A)nS,t11 (0  ,r-l;k-l,2:A) 
n.r  n,r+l 


=  nS,t(0,r-l;k-l:A)nS,t1 . (0 ,r-l ;k, 2 : A) 
n , r  n , r+1 

where 


(8.5)  nS,tL1  (0,r-l;k-l;k+l:A)  = 
n,r+l  . 


P 


r-t+1, . . .  ,n-s-t  ,n-s+l, . .  n 

0  , . . . ,n-k-2 ,n-k ,n-k+2 , . . . ,n-r+l 


Proof . 


Consider  the  following  matrix  of  order  (n-r+1)  x  (n-r+2) 


32. 


(8.6) 


1  ^  (1-X)  0  0 


1  /t-t+1. 

1  (  i  ) . (1-A)  0 


,n-s-t 


.n-s+1. 
^  1  ' 


(^f1) 


1  (") 

M  II 
f(D  f (2) 


..  ( 


0 


,0 


1  (  !  > . (1-A)  0..0 _ 0, 


,0 


(1-A)  0 - 0 


(1-A) 


n  )(  n  )(n')('  n  )  (  n  ^ 
n-k-2;  n-k-l' n-k' ^n-k+l'  ^n-k+2' 


(  n  ) 

*  Vr+l' 


(n-k-1) 

.  .  r  a 


II 

(n-k) 


•  f 


(n-r-1) 


Let  f^  denote  the  i-th  column  (i  =  1 , 2 , . . .  , n-k-1)  ,  a,  b,  c 

(  Q\ 

denote  the  (n-k),  (n-k+1)  and  (n-k+2)-th  columns  respectively  and  r 

denote  the  v-th  column  (v  =  n-k+3 , . . . ,n-r+2)  of  (8.6)  as  shown  above. 

T 

If  we  set  d  =  (0,...,0,1)  we  obtain  (8.3)  from  (6.2),  and  if  we  take 
d  =  (1,0,...,0)T  we  get  (8.4).  □ 

Proof  of  theorem  8.4  The  proof  follows  by  induction  on  k.  Observe 

that  since  n S ’ t (0 ,r-l;r-l :A )  =  nS,t(A),  theorems  8.2  and  8.3  assert  that 
n,r  ’  n,r 

theorem  8.4  is  true  for  k  =  r-1  and  r. 


33. 


Let  us  suppose  that  the  assertion  is  true  for 
k  =  r-1, ...,£.  W.  1.  o.  g.  we  may  assume  that  r  is  even,  and  let  us 


denote  the  zeros  of  II® ’^(0,r-l;£:A)  by  {y^}  ^  so  that 


n,r 


(8.7) 


UK  (£) 


(£-lK  (£) 


Y  1  "  Y'l  "'<  r2'<-<  Yi-7<  Y'H r+L  0  " 


(£)  (Jl)  (£-1) 

<v  <  <vx/<vv  ' 

Y  £-r+2  Y  d  Y  d 


It  follows  from  (8.3)  and  (8.4)  that 


(8.8)  n®’ho,r-l;  £+1:  Y  (f  )n®  },ho,r-2;H-2.2:Y<f  ) 

II  y  -I  J-  II  I  j  I  1 


=  -ns,t(0,r-l;£-l:y(£))ns  ] ’ C (0 ,r-2 ; £-1 ,2 :  y  (£) ) 
nrr  i  n-l,r  i 


and 

(8.9) 


nf ’I(0,r-l;S,+l;Y(^))n®’^+1(0,r-l;J,-l,2:Y^h 


n,r 


rs,t 


(£).TTs,t 


(£) 


-II_  ’^(0  ,r-l;£-l;y  .  )n  ’  (0,r-l;£,2:y  .  ) 

1  H  y  IT  _L  i 


n,r 


(i  =  1,2,. . .  ,d  ,  r  <_  £  <_  j-2)  . 


rS,t 


Since  11“  *  w  n  (0  ,r-l ;  £-1, 2  :y  ^)  and  ns,t  (0,r-l;£,2:y^) 
n,r+l  ’  l  n,r+l  i 


are 


q— 1  i~  . . —  ■  ■  ■  ...  (  Q  \ 

positive  for  i  =  £-r+2,...,d,  and  ^ (0 ,r-2 ; £-2 , 2 :y  _/  )  and 


.(£) 


IIs  |,t:(0,r-2;£-l,2:y  are  positive  for  i  =  1, 2 , . . . , £-r+l ,  it  follows 

n  ”*  j  it  i 

that 


(8.10)  Sgn  n®,^(0,r-l;£+l:yv^)  =  -Sgn  n^’^(0,r-l;£-l:y ) 


(£) 


rs»t 


(£) 


(i  =  1,2,. ..,d). 


TS. 


rs,t 


Hence  II  ’  (0,r-l; £+l:X)  has  exactly  one  zero  in  each  of  the  intervals 


(&)  on 

(Y  i  »  Y  )  (i  =  1, 2, . . . ,d-l) .  Clearly  (8.7)  shows  that 


rs,t 


II  *  (0,r-l;£-l:Y 


n,r 


00 

£-r+l 

£-r+l 


)  >  0,  since  IIs  ’ (0  ,r-l ; £-1 : 0)  >  0.  Hence 


n,r 


IIs  ’  (0  ,r-l ;  £+1:  y  ^  )  <  0.  But  then  IIs  ’  ^  (0  ,  r-1 ;  £+1  :  0)  >0.  Hence 

n ,  r  ■**  ^  v- 


rs,t 


n,r 


00 


.00 


the  zero  of  II  *  (0  ,r-l ; £+1 :  X)  between  (y  ,  y  9)  actually  lies 

n  ^  it  x.—V’t  xz-rT  z 

(£)  s  t  - 

in  (y 0).  Since  II ^  (0  ,r-l ;  £+1:  X)  is  a  reciprocal  polynomial  it 


n,r 


.(£) 


must  also  have  a  zero  in  (-°°,Y  ^  ).  The  assertion  follows  by  induction, 


An  analogous  method  gives  the  following  generalisation  of 
theorem  7.5. 


Theorem  8.6  Let  r  >  2,  and  let  n  =  2m-l  be  an  odd  integer  such  that 
m  >  r+1.  Then  for  each  integer  k  with  r-1  <  k  <  n-r-1, 

nn  r,^(0,r-l;k:X)  has  real  simple  zeros,  (k-r+2)  of  which  are  of  sign 
n,r 

r — l  r 

(-1) _ and  (d-k+r-2)  of  sign  (-1)  .  Furthermore,  the  zeros  of 

j[n  r»^(o  r-l;k:X)  and  IIn  r,^(0,r-l;k+l:X)  are  separated  for 
n,r  n,r 

k  =  r-l,r , . . . ,n-r-2. 


9.  Discussion  of  the  C.L.I.P. 


From  the  results  of  §8  we  see  that  if  n  =  2m-l  is  odd,  the 


s  t 

zeros  of  the  polynomial  II  *  (0,p;k,q:X)  are  real,  simple  and  not  equal 

n ,  r 

to  ±1  for  each  of  the  following  cases: 


(9.1)  t  =  0,  k  =  n-r-q+1,  n  >  2r+l, 

(9.2)  t  and  q  are  both  even  and  0  <_  t  <  r  <  k+q  s  <_  n-r+1, 


. 


(9.3) 


t  even,  q  =  1  and  0  <_  t  <  r-1  <  k  s-1  <_  n-r  with  (k-r)  odd, 

(9.4)  t  =  1,  q  =  1,  s  =  n-r  and  0  <  r-1  <  k  <_  n-r-1  with  (k-r)  even. 

The  method  of  Chapter  I  leads  to  the  following  theorem. 

Theorem  9 . 1  Let  p,q  satisfy  (1.1)  and  suppose  that  p,  q,  k,  s,  t 
satisfy  one  of  the  conditions  (9.1)- (9. 4) .  Then  given  r  bi-infinite 

sequences  (1.5)  satisfying  (1.6)  the  C.L.I.P,  of  type  (0,p;k,q)  has  a 

s  t 

unique  solution  S(x)  e  -^m-l  r  suc^  t^iat 

(9.5)  S (x)  =  0 ( | x | Y)  (p  e  0 ,p;k,q) 
as  jxl  00 . 

If  one  of  the  zeros  of  the  characteristic  polynomial  is  1 
(or  -1)  the  corresponding  C.L.I.P.  does  not  have  a  unique  solution. 
Indeed  if  S^(x)  is  an  E.S.  corresponding  to  the  E.V.  1  (or  -1),  then 
S^(x)  +  S(x)  is  a  solution  of  the  C.L.I.P.  satisfying  9.5,  if  S(x)  is 
a  solution  of  the  same  problem  satisfying  (9.5). 

Remark .  We  are  concerned  here  only  with  a  particular  class  of  C.L.I.P., 

- -  s  t 

namely  problems  of  type  (0,p;k,q)  in  S  *  in  which  only  one  gap  is 

n ,  r 

allowed  in  the  interpolating  conditions.  Even  in  this  particular  class 
we  do  not  completely  consider  all  the  cases.  Conspicuously  absent  is 
the  case  when  q  and  t  are  both  odd.  The  only  problem  of  this  type 
which  we  have  considered  is  given  by  condition  (9.4). 


CHAPTER  III 


HANKEL  DETERMINANTS  OF  EULER-FROBENIU S  POLYNOMIALS 

S  t  1  "  ■  ■■  ■  ■ 

The  characteristic  polynomials  IIn * r (0 , p ;  k,q:A)  play  a  funda¬ 
mental  role  in  the  corresponding  C.L.I.P.  Since  much  is  known  about 

the  Euler-Frobenius  polynomial  II  (A)  it  is  natural  to  search  for 

n 

s  t  - - 

relations  between  II  ’  (0,p:k,q:A)  and  II  (A).  It  turns  out  that 

n,r  r  n 

II  (A)  is  related  to  a  Hankel  determinant  of  II  (A)  .  We  shall  also 
n,r  n 

state  without  proof  an  expression  for  the  Hankel  determinant  of  expo¬ 
nential  Euler  polynomials.  This  relation  also  occurs  in  an  interpo¬ 
lation  problem  with  cardinal  splines,  but  we  shall  not  dwell  upon  it 
here. 


10.  A  relation  between  II  (A)  and  II  (A). 

n  ,r  n 


Let  us  denote  by  H  , (a  )  the  Hankel  determinant  of  order 

J  r+1  n 

(r+1)  x  (r+1)  given  by 


(10.1)  Hr+l(an) 


a  i 

n 

n-1 

n-r 

n-1 

a  „ 
n-2 

n-r-1 

a  a  ,  .  a  0 

n-r  n-r-1  n-zr 


Determinants  of  this  type  in  which  the  entries  a^  are  orthogonal  poly 

nomials  have  been  studied  by  several  mathematicians  (see  [1]  and  [3]) 

If  a  =  II  (A)/n!  we  have  the  following  relations, 
n  n 


36 


Theorem  10.1.  Let  n  and  r  be  non-negative  integers  such  that 


n  ^  2r+l.  Then 


n  (A) 

(10-2)  Hr+1(  ^T-> 


^I±i] 

(-1)  2  2  !3! . . .r! (1-A)r(n_r)n  (A) 

_ n , r+1 

n! (n-1) ! . . . (n-r) ! 


The  proof  of  theorem  10.1  depends  on  the  following  lemma,  which  is  a 
particular  case  of  a  more  general  identity  on  symmetric  determinants 
(see  [9],  p.  372). 


Lemma  10.2.  For  n  >  2r+l  we  have 


(10.3) 


[Hr  <Vl>]  -  Hr-l(an-2>Hr+l<0 


n 


Proof  of  Lemma.  The  lemma  is  easily  established  using  the  determin- 
antal  identity  (6.2)  on  the  following  vectors:  c  and  d  are  the  first 
and  last  columns  of  (10.1)  respectively,  f^  the  (v+l)-th  column 
(v  =  l,2,...,r-l)  and  a  =  (1 ,0  , . . .  ,0) T ,  b  =  (0,0, . . . ,0,1)T. 

Proof  of  Theorem  10.1.  The  proof  will  be  carried  out  by  induction  on 
r.  First  of  all,  let  us  observe  that 

n  (A)  n  (A) 

H  ( - - — )  =  - —  for  all  n  >  1 

1  n!  n!  — 


and  by  (7.1)  of  lemma  7.1,  with  r  =  2,  we  have 

n  (X)  (-i)(i-x)n_1n  (A) 

H2(-^r)  =  - n !  (n-1)  !  ■’ - -  for  311  n  -  3  • 

Let  us  suppose  that  (10.2)  holds  for  H^(II  (A)  /n! )  (k  =  0,1,..., r)  and 
we  shall  prove  that  (10.2)  also  holds  for  H  ^(n  (A) /n! ) . 

Using  (10.3)  after  some  calculations,  we  have,  for  n  >_  2r+l, 


. 


Vx)  [ir] 

(10-4)  V2,r-l(X) 


2 !  3 !  .  .  . (r-2) ! (r-1) ! (1-A) 
n!  (n-1)  !  .  . .  (n-r)  ! 


r(n-r) 


{n(n  (X)) 

n-l,r 


(n-r)  II  (A)  (X)}  , 

n ,  r  n-  Z ,  r 


from  which  we  obtain  (10.2),  with  the  help  of  (7.1). 


g  £ - - 

11.  Some  relations  between  II  ’  (0,p;k,q;A)  and  II  (X)  . 

n.r  n 


Relations  similar  to  (10.2)  involving  other  polynomials 

g  £  — ■ ..  — — — 

II  ’  (0,p;k,q:X)  can  be  derived  by  the  same  technique.  Set 
n  >  r 


a 

n 

a  n 
n-1 

•  •  • 

an-r+l 

a  i 

n-r-k 

a  i 
n-1 

a  0 
n-2 

•  •  • 

a 

n-r 

an-r-k-l 

(11.1)  Hr+1(0,r;k:an)  = 


and 


(11.2) 


H?+l(an) 


a 

n-r 

a  i 

n-r-1 

’ *  *  an-2r+l 

an-2r-k 

a 

n 

a  i 

n-1 

‘ * '  an-r+l 

a  i 

n-r-1 

a  , 

a  o 

•  •  •  3 

a  „ 

n-1 

n-2 

n-r 

n-r-2 

3  3  •••  3  ~  3  ^ 

n~r+l  n-r  n-2r+2  n-2r 

n-r-1  n-r-2  n-2r  n-2r-2 


then  using  a  method  similar  to  that  of  Theorem  10.1  we  have 
Theorem  11.1.  Let  n,r  be  positive  integers.  Then 


(11.3) 


n  (x) 


n 


Hr+i(°>r;k:"~T'  ■) 


n! 


[I±l] 

(-1)  2  1! 2! . . .r(l-l)r(n~r)~k 


n 


n! (n-1) ! . . . (n-r+1) ! (n-r-k) !  n,r+l 


(0,r;r+k: A) 


(n  >_  2r+k+l,  k  >_  0) 


n  (X) 


[I±I] 

(-1)  2  l!2!...(r-l)!(r+l)!(l-l)r(n~r)~2 

n! (n-1) ! . . . (n-r+1) ! (n-r-1) ! 


I8 

n,r+l 


X  n&  ...  (A)  (n  >  2r+3)  . 


p 

where  the  polynomials  II  ^  r+^(X)  are  given  by  (3.8)  and  the  polynomials 
^n  r+p(0 >r »r+k:X)  are  given  by  (5.3)  and  (5.3a). 


12.  Hankel  determinant  of  exponential  Euler  polynomials. 

Using  the  same  technique,  a  more  general  form  of  the  relation 
(10.2)  can  be  obtained  in  another  direction.  More  precisely,  if  A  ^(x;A) 
are  the  exponential  Euler  polynomials  (see  §13),  then  for  r  =  0,1,2,..., 
n  >_  2r+l, 


(12.1) 


(x;  X) 
n ! 


) 


[__]+(r+i)n 

(-1)  _ 1!  2!  — r! 

n ! (n-1) ! . . . (n-r) ! 


Vr+l(x;X)! 

(i-x)n~r  j 


9 


where 


40. 


1 

1 


(j) . (r-i)(1_x)  o  o 


<rjb 


(r*b 


(l-X)O 


0 

0 


(12.2)  A  (x;  A) 
n ,  r 


1  (*) . 

n-r+1  .n-r+1.  n-r 
x  (  1  )x  .  1 

n-r+2,n-r+2.  n-r+1 
x  (  1  )x  . 


•  •  Ci>  <llx> 
0 . 0 

1  0 . 0 


n-1 

x 


n 


Ci)x  1 


In  view  of  the  relation  A  (0 ;  A)  =  H  (X)/(A-l)n  (see  §13),  it  is 

n  n 

clear  that  (10.2)  is  a  particular  case  of  (12.1)  when  x  =  0.  A  full 
discussion  of  the  relation  (12.1)  and  the  corresponding  interpolation 
problems  will  be  discussed  elsewhere  in  a  joint  paper  with  A.  Sharma 


and  J.  Tzimbalario. 


. 


CHAPTER  IV 


EXPONENTIAL  HERMITE  EULER  SPLINES 


In  a  very  interesting  paper  [14]  Schoenberg  studied  the  Cardinal 
splines  that  interpolate  the  functions  X  at  the  integers,  where  X  is  a 
complex  number  (see  also  [15]).  He  called  them  the  exponential  Euler 
splines,  because  they  are  "periodic  extensions"  of  the  exponential  Euler 
polynomials  A  (x;X)  generated  by  the  relations 


n 


X— 1  xz 

-  e 


00  A  (x;X) 


X-e 


V  n  n 

=  l  — sr  2  • 

n=0 


This  chapter  deals  with  Cardinal  splines  that  interpolate  the  function  X 

as  well  as  its  consecutive  derivatives.  The  problem  is  to  find  a  spline 

S  (x)  e  S  such  that  S^P\v)  =  (log  X)P  XV  for  all  integers  v 
n,r  n,r  n,r 

(p  =  0,1, ... ,r-l) . 


13.  The  polynomial  A  (x;X)  . 

t ,  s 

Let  n,r  be  positive  integers  such  that  n  >_  2r-l  and  let 
s  =  0,1,..., r-1  be  a  fixed  integer.  Let  us  define  the  polynomials 
A  (x;X)  for  0  <  x  <  1  as  follows: 

n  jT  jS 


-  41 


. 

_ 


42. 


(13.1)  (-l)faA. 


n,r,s 


(x;X)  = 


An(x;  A) 


n! 


A  . (x; A) 
n-1 


(n-1) 


An(0;A) 

n! 


A  ..  (0;  A) 
n-1 _ 

(n-1) ! 


Vr+l(x;X)  Vr+l(0;X) 


(n-r+1) ! 


(n-r+1)  ! 


Vs+l(0»X> 

(n-s+1) ! 


An_s(°;X) 


(n-s) ! 


A  0 (0 ; A) 

n-r-s+2 _ 

(n-r-s+2) ! 


Vs-l(0;X)  Vr+d0^ 


(n-s-1)  ! 


A  0 (0 ; A) 
n-s-2 

(n-s-2)  ! 


A  (0; A) 
n-r-s 

(n-r-s) ! 


(n-r+1) ! 


(n-r) 


An-2r+2(Q;X') 

'  (n-2r+2) ! 


where  the  exponential  Euler  polynomials  A^(x;A)  = 

xn  +  (^l)a1(A)x11  1  +  ...  +  a  (A),  and  A  (0;A)/n!  =  a  (A)/n!  = 
ly  1  n  n  n 


II  (A)/n!(A-l)n  (see  [14]).  It  follows  from  (10.2)  and  (13.1)  that 
n 

/  A  (0; A) 


(s) 


n,r,s 


(0; A)  =  H  (- 


n 


n! 


-)  = 


(13.2)  < 


[f  ]+(r-l) (n-r+1) 

(-1) 


1!  2! . (r-1)  !  n  (A) 


na.r 


n!  (n-1)  ! . (n-r+1)  !  (A-l) 


n-r+1 


Using  the  relations  A  (x;A)/n!  =  A  (x; A) / (n-1) ! ,  and 

°  n  n-i 

A^P^  ( 1 ; A)  =  AA^(0;A)  (p  =  0,1, . . .  ,n-l)  ,  it  is  easy  to  check  that 
n  n 

A  (x; A)  satisfy  the  relations 

n,r,s 


43. 


(13.3)  A(p)„  Jl;X)  =  XA(p)  (0;X)  (p  =  0 ,1, . . .  ,n-r)  , 

n,r  ,s 


n,r,s 


^  A(p)  ( 1 ; X)  =  A(p) 


n,r,s 


n,r  ,s 


(0 ; A)  =  0  (p  =  0,1, . . . ,r-l,p  4  s)  ,  and 


(13.4)  / 


(s) 

1  n,r , s 


A  (  0 ;  X  ) 

(0;X)/H.  -S-; - 1  =  1 


n! 


provided  X  ^  1  and  X  is  not  a  zero  of  II_  ^(X),  an  assumption  which  we 
shall  impose  throughout  this  chapter. 


n,r 


Observe  that  when  r  =  1  (in  which  case  s  can  take  only  the 

value  zero),  A  ,  (x;X)  =  A  (x;X) . 

n ,  1 , 0  n 


14.  The  spline  S  (x;X). 


n,  r ,  s 


Let  us  define  a  function  S  (x;X)  such  that 


n,r  ,s 


(14. 1-)  S  (x;X)  =  A  (x;X)/H 
n,r  ,s  n  j  r  j  s  ic 


\(0;X)' 


n! 


(0  <  x  <  1)  ,  and 


S  (x+l;X)  =  XS  (x;X)  for  all  real  x. 

n,r,s  n,r,s 

It  follows  from  (13.3)  and  (13.4)  that  S  (x;X)  e  C  (-°°,co)  ,  and 

n ,  r ,  s 

/ 

S ^p)  (v:X)  =  0  (p  =  0,1, . .  .  ,r-l,p  4  s)  and 

n,r,s 


(14.2)  l 


S(s)  (v:X)  =  XV  (v  =  0,±1,±2,. . .) . 


n,r  ,s 


so  that  it  is  a  Cardinal  spline  with  integer  knots  of  multiplicity  r, 
and  belonging  to  the  class  =  (S(x)eS  :  S^(v)  =  0  V  integer  v, 

I"l  j  XT  XX  y  L 

p  =  0,1,... ,r-l;  p  4  s>  . 


Now  let  us  set 


r-1 

(14.3)  S  (x,*)  =  I  (log  *)S  S  (x;^)  (x  e  R) . 

s=0  n,r,S 


n,r 


The  following  theorem  is  an  easy  consequence  of  (14.2). 


44. 


Theorem  14.1.  The  splines  ^(x;A)£S^  ^  and  satisfy  the  following 


intepolatory  conditions: 

(14.4)  S^P^(v;A)  =  (log  A)p  AV  (p  =  0,l,...,r-l,  and  V  integers  v) . 

n  y  r 

15.  Convergence  of  exponential  Hermite-Euler  splines  when  r  =  2  ♦ 

When  r  =  1,  S  ..  (x;A)  E  S  (x;A)  are  the  exponential  Euler 

n ,  I  n 

splines  considered  by  Schoenberg  [14]  who  proved  that  lim  S  (x;A)  =  A 

n 

n-*® 

uniformly  for  all  x  belonging  to  a  finite  interval,  when  A  is  a  non¬ 
negative  complex  number.  When  r  =  2  we  have  the  following 


Theorem  15.1.  If  A  is  not  a  positive  complex  number,  then 

(15.1)  lim  S (p^  (x; A)  =  (log  A)P  AX  (p  =  0,1) 

n,  z 

n-*» 

uniformly  for  x  belonging  to  a  finite  interval. 

We  shall  show  that  (15.1)  is  an  easy  consequence  of  the 
corresponding  results  on  S  (x;A).  Let  us  define  two  functions 

in.  y  Z_  y  S 

,,  -  ,  s  ,  X  -TTix  sin  TTX 

(15.2)  a(x)  =  A  e  -  , 

TT 

(15.3)  g(x)  =  AX  -  (log  A)  a(x) 

Theorem  15.1  clearly  follows  from  the  following 

Theorem  15.2  .  Let  A  =  |  A  [  ela  .  If  0  <  a  <  2tt  ,  the  following  relations 
hold  uniformly  for  all  x  belonging  to  a  finite  interval: 

(15.4)  lim  S(p^  n(x;A)  =  g(p)(x)  (p  =  0,1)  , 


and 

(15.5)  lim  S(pJ  ,(x;A)  =  a(p)(x) 

n,  2 , 1 


(p  =  0,1)  . 
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ia 


Lemma  15.3.  Let  X  =  |X|  e  and  X  =  log  |X|  +  i(a+2iTk)  .  The  following 

1\ 

relations  hold  uniformly  for  all  x  in  [0,1]: 


(15.6) 

lim 

n+1 

A 

n-*» 

U 

(15.7) 

lim 

n+1 

A1 

n-x° 

(15.8)  lira  X 

n-xc 


A  (x;X) 
n _ 

n! 


=  (X-l)X  1XX  (-7T  <  a  <  it  )  , 


A  (x; X)  X  A  (x;X) 
n-1  on 


(n-1) : 

(-it  <  a  <  0) 


n! 


}  =  (X-l)A_1AXe^7rix  (2Tri) 


n+1 

-1 


{ 


A  1 (x;X)  X  A  (x;X) 
n-i  o  n 


(n-1) ! 

(0  <  a  <  tt)  . 


n ! 


=  (X-l)X  1XXe  2^1X(-27ri) 


Proof .  Using  the  expansion  (see  [14]  p.  399) 


(15.9) 


A  (x; X) 


n 


n! 


=  (X-l)X  1Xx 


k=-° 


2TTkix  ,,  n+1 
e  /xk 


(0  <  x  <  1)  , 


we  have 
(15.10)  X 


n+1  An^x>x)  -l  x  r  2-rrkix  .  ,  . 

=  (X-l)X  X  l  e  ( X  / X^) 

k 


n+1 


o 


n. 


Since  |X  |  <  |x^|  V  k  ^  0  ,  (15.6)  follows  from  (15.10). 

Also  from  (15.9)  we  have 

A  .  (x;  X)  X  A  (x;  X) 

<«-u>  xf  {-^rbryi - ’Hh — »  ■  (x-1)x  x  * 

v  /  ,  w\n+1  2ukix 

x  l  e 

k^0  k 

If  -tt  <  a  <  0,  |  X  |  <  IXJ  V  k  ^  0,1,  and  (15.7)  follows  from 

(15.11)  .  The  limit  (15.8)  is  proved  in  the  same  way.  □ 


Proof  of  Theorem  15.2.  We  shall  prove  only  the  relation 


(15.12) 


lim 

n-*» 


S  0  n (x; X)  =  3(x)  . 
n,  2 ,0 


The  rest  are  proved  in  the  same  way. 


j iSL '  •  <  >”  +  M  3‘ 1  *  / 


■ 


If  0  < 


We  can  write  An+1ArVj~1A  0  ~(x;A)  = 

o  -1  n,2,0 


n+1  An(x;A) 


o  n! 


A 


n+1  An-l(0;X) 
o  (n-1) ! 


A  (x;A) 
,  n+1 r  n-1 _ 

-1  1  (n-1)! 


A  A  (x; A)  A  _(0;A) 

on  ■)  , n+1  r  n-2 

}  -1  1  (n-2)! 


n! 


A  A  ,  (0 ; A ) 
o  n-1 _ 

(n-1)  ! 


a  <_  ir,  it  then  follows  from  (15.6)  and  (15.8)  that 


(15.13)  An+1Ant1A  0  n(x;A)  (A-1)2A  2AX(A  -A  )  (A  -A  e  2TTlX) 
o  -1  n,  2  ,U  -1  o  -1  o 


Hence 


x  -2TrixN 

(  x  (X_i"X0e  > 

lim  Sn  ?  n(x;X)  =  — 7\ - T~\ - 

n  >  2 ,0  (A_1  -Aq) 


(15. 14)  \ 


n-x» 


=  Xx{1  _(log  AX!~e  ,.2TriX) 

2Tri 


}  (0  <  a  <  it)  . 


( 


Similarly 

,.  „  (  ,  xf  (log  »(l-e27T1X) , 

llm  Sn,2,0(x;X)  =  X  {1  +  - 2ifi -  1  " 

n->-oo  ’  * 


(15.15) { 


“  2  TT 1  X 

,x  2tt1x r .  (log  X  +  2iri)(l  -  e  x  ) 
A  e  {1 - - 2iri  ' 


}  ( — TT  <  a  <  0) 


Combining  (15.14)  and  (15.15)  we  obtain 
(15.16)  lim  S 


n-H*> 


n,2,0 


(x; X)  =  XxU  -(loS  «l-e': 2*1X) 


} 


when  A  =  |A|eia  for  0  <  a  <  27T,  from  which  (15.12)  follows. 


□ 


Remarks.  For  the  convergence  of  the  exponential  Hermite-Euler  splines 


we  deal  only  with  the  case  r  =  2.  The  general  results  can  perhaps  be 
obtained  by  the  same  method. 
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16.  Exponential  Euler  g-splines. 

The  techniques  of  the  preceding  sections  can  be  applied 

without  difficulty  to  study  C.L.I.P.  for  the  function  A  .  For 

instance,  let  us  consider  the  problem  of  finding  a  spline 

S8  0(x)  e  S8  such  that  S8^  (v)  =  (log  A)P  AV  (p  =  0,2,  and  V 
n ,  Z  n,  Z  n,Z 

integers  v) . 

In  this  case  we  introduce  the  polynomials 


(16.1) 


,8  o 

n ,  2 ,0 


A°  „  „ (x; A) 


A  (x; A) 
n _ 

n! 


A  0 (x; A) 
n  z 


(n-2) ! 


|An(x;A) 


,8 
n ,  2 , 2 


A°  «  o(x;A) 


n! 


A  9 (x; A) 
n-Z 


(n-2) ! 


A  9 (0 ; A) 
n-2 _ 

(n-2) ! 


A  ,  (0 ; A) 
n-4 

(n-4) ! 


An(0;A) 

n! 


An-2(0;X) 

(n-2) ! 


and  define  the  splines 


(16.2)  S8  _  (x; A)  =  A8  (x;A)/H8 

n, z, s  n , z , s  z 


/a  (0;A)\ 
n 


\ 


n! 


for  x  e  [0,1],  and  S8  0  (x+l;A)  =  AS8  (x;A)  for  all  x  e  R  (s=0,2) 

n,z,s  n , z , s 


Then  S8  _  (x;A)  e  S8  (s  =  0,2),  and  we  have 


n,2  ,s 


n,r 


Theorem  16.1.  The  spline  functions 


(16.3)  S8  „ (x; A)  =  S8  ft(x;A)  +  (log  A) 2  S8  2(x;A) 
n,Z  n,z,u  ii ,  i.  y  & 


belong  to 


and  satisfy  the  following  interpolating  conditions 


(16.4) 


sg(p> 

n,2 


(v ; A)  =  (log  A) P 


(p  =  0,2,  and  V  integers  v) . 


CHAPTER  V 


FOURIER  TRANSFORMS  OF  B-SPLINES  AND  FUNDAMENTAL  SPLINES 

Fourier  transforms  of  B-splines  and  fundamental  splines  provide 
a  powerful  tool  in  the  C.I.P.  (see  [7],  [10],  [12],  [15]).  As  an 
application  of  the  exponential  Hermite-Euler  splines,  we  shall  compute 
the  Fourier  transforms  of  B-splines  and  fundamental  splines  belonging 
to  S  _  for  (0,1)  interpolation.  The  same  method  can  be  applied,  with 

XT  y  Z 

suitable  modifications  to  the  general  C.L.I.P. 

17.  Fourier  transforms  of  L.  ..  _  (x) . 

 zm-l , Z ,  s 


Let  us  recall  that  if  a  function  f (x)  satisfies  the  condition 


(17.1)  f^(x)  =  0(|x|Y)(p  =  0,1)  as  |  x|  ->  00  ,  y  >  0  , 


then  there  exists  a  unique  spline  ^2m_-^(x)  £  5 


such  that 


f(p)(v)  (p  =  0,1)  V  integers  v 


and  S 


2m-l 


(x)  is  given  by  the  Hermite  interpolation  formula 


00 


00 


(17.3) 


f  (v) 


where  L 


2m- 1 , 2 , s 


(x)  (s  =  0,1)  is  uniquely  determined  by  the  conditions 


(17.4)  Lkp;  _  (v)  =  0  V  integers  v  ±  0  (p  =  0,1)  , 

2m-l,2,s 
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Theorem  17.1.  The  Fourier  integral  representations  of  the  fundamental 


splines  are  given  by 


(17.6) 


__L 

2tt 


L2m-l,2,0(x) 


\, 2  sin  U/2v^m  ,  N  ,2  sin  U/2N^m  ^  /  s 

iux[( - u - ■>  p2m-2(u)-( - U - }  p2m-l(u) 


du. 


H2(p2m(u)) 


L2m-l,2,l(x) 


(17.7)  < 


1 

2tt 


u  2m 


f/2  sin  /2*  ,  x  f 

luxl  u  2m-l  - 


.  u.  2m- 1 
2  sm  /2n  ,  N 

- )  P,  (u) 

u  2m 


du, 


(2  sin  U/2)H2(p2m(u)) 


V 


where  (using  the  notation  of  Schoenberg  [12],  p.  178) 

00 

(17.8)  p  (u)  =  (2  sin  U/2)n  V  l/(u+2Trj)n  . 
n 

J=-°° 

Proof.  Let  S„  ..  (x;e1U)  (s  =  0,1)  be  the  exponential  Hermite-Euler 

-  2m-l,s 

splines  of  Chapter  IV  with  X  =  e1U(0  <  u  <  2tt)  .  Following  Schoenberg 
[16],  it  is  easy  to  show  that  the  functions 


(17.9) 


_1_ 

2tt 


2tt 


S9  ,  9  (x;e1U)du  (s 
2m- 1 , 2 , s 


0 


0,1) 


are  spline  functions  belonging  to  2*  Furthermore  the  functions 

(17.9)  satisfy  conditions  (17.4)  and  (17.5)  in  view  of  the  inter¬ 
polating  conditions  S^2^_^  2  s(v;e^U)  =  e^  U  (s  =  0,1,  v  =  0,±1,±2, . . .) , 

and  S(p^  .  0  (v;e1U)  =  0  (p  ^  s,  v  =  0 ,±1,±2 , . . . ) .  Hence  from 

2m-l,2,s 

uniqueness 


51. 


(17.10)  L  .  (x) 

2m-l,2,s 


Jl 

2tt 


2ir 


S2m-l,2,s^X’e  ^  du  °,1') 


0 


From  (14.1)  we  have  for  0  <  x  <  1 


(17.11) 


S2m-l,2,0('X;e  ^ 


iux  A2m-1 , 2 ,0 (x>  elU) 


A2m-l,2,0(‘0;e  ) 


Using  (13.1)  and  (15.9)  we  can  write 


(17.12) 


e"1UXS 


2m 


,  iu,  A2m-l(x;u) 
-l,2,0(x’e  )  &2m_1(0;u) 


where 


(17.13)  ^2m-l^x’u^  = 


I  e2”kiX/(u+2,k)2m  l  1/  (u+2iTj )  2m-2 


k=-° 


l  e2lrkix/  (u+27Tk)  zm-±  l  l/^Trj) 
k=-°°  j=-3° 


2m- 1 


. x 2m- 2 


Clearly  A„  .  (x;u)  is  a  2ir-periodic  function  of  u.  From  (17.10),  (17.12) 
2m- 1 

and  (17.13),  on  interchanging  the  order  of  integration  and  summation 
(can  be  justified  by  Fubini  theorem)  we  obtain 


(17.14) 


—  I 

2tt  l 


k=-c 


L2m-l,2,0(x) 


r27T  i(u+27rk)x 
e  _ 

0  A2m-l(0;u) 


2m 


.  N  2m-l 


1/ (u+27Tk)  \  1  /  ( U+  2 7T  j ) 

-1  =—00 


l/(u+27Tk)2m  1  £  l/(u+2TTj)2m  2 

-i=_oo 


du. 


Since  A  (0;u)  is  a  2-rr-periodic  function  of  u,  by  a  change  of  variable, 
2m- 1 


we  obtain 


52. 


(17.15) 


L2m-l,2,0(x) 


2tt 


1UX 


l 


2m  .L  (  0  ..2m- 2 

L  u  1=-°°  (u+2Tr-|) 


00 

y 

2m-l  . ^ 


A2m-l(0;u) 


U_ j=-oo  (u+2TTi) 


2m- 1 


du 


Mutatis  mutandis,  we  have 


(17.16) 


L2m-l,2,l(x) 


l 

2tt 


1UX 


I 


2m  #  /  i  a  •  \  2m  1 

u  j=-°°  (u+2ttj) 


00 

—  y. 

m-l  .  ^ 


2m- 1  . L  .  ..2m 

u  j=-°°  (u+2ttj) 


du 


A2m-l(0;u) 


Straightforward  calculation  using  (17.8)  shows  that 


(17.17)  A  ,(0;u)  =  (2  sinU/2)  4m+2 
2m- 1 


p2m(u)  P2m-2(U)  '  pL-l(u)] 


=  (2  sin  U/2)  4m+"2  H2(p2m(u))  , 


and  (17.6),  (17.7)  follow  from  (17.15)  and  (17.16)  respectively. 


□ 


18.  Fourier  Transforms  of  B-splines. 


The  B-splines  for  C.H.I.P.  have  been  studied  by  I.  J.  Schoenberg 
and  A.  Sharma  [17].  For  the  case  when  r  =  2  these  B-splines  N  (x) 

o 

(s  =  0,1)  are  defined  by 

(m-2) 


(18.1)  N  (x) 
s 


I 


v=-(m-2) 


Cv  ^2m-l,2 


S(X_V)> 


53. 


where  c  are  the  coefficients  of  the  Euler-Frobenius  polynomials 


v 


n2m-l,2(X)  given 


(18‘ 2)  n2m-l,2(X)  C-(m-2)  +  °-(m-3)X  +  *  *  *  +  C0X  +  •••  + 

c t  9a2m-4  (c  >  0)  . 

(m-2)  o 

For  the  B-splines  N  (x)  (s  =  0,1)  we  have  the  following  representation 

s 

theorem. 

Theorem  18.1.  The  Fourier  integral  representations  of  the  B-splines 
are  given  by 


m 


(18.3) 


H  00  -  X 


1UX 


(2  sin  U/2) 2  ,  s  (2  sin  U/2) . 

- 2S -  P2m-2(U)  '  "  2m'-l  P2m-l(u) 


u 


u 


du. 


m 


(18.4) 


N1(x) 


(2m-l)  !  (2m-2)  !  (-1)  i 

2tt 


rco 


1UX 


(2  sin  u/2) 

2m 

u 


p2n,-l(u) 


u 


2m- 1  2m 


PoJU) 


du. 


Proof.  Using  (17.15)  and  (18.1)  we  have 


,00 


(18.5)  N  (*>  -  ^ 


xux 


I  (rn-2) 

l 

^-(m-2) 


c  e 
v 


— xvu 


_A_  y 1 — 

u2m  -°°  (u+2tt  j )  2m 2  u2m  J'  (u+27tj) 


00 

—  y  __i — 

i-l  L  ~  .  x2m-l 


du  , 


A2m-l(0;u) 


54. 


whence 


(18.6)  N  (x) 
o 


_1 

2tt 


iux  -(m-2)  iu  ,  iu. 

6  6  n2m-l,2(e  5  X 


CO 

—  y  — 

L 


oo 

—  y 

>m—  1  ^ 


2m  L  (  ..2m- 2  2m-l 

u  -°°  (u+2ttj  )  u  — °° 


,  10  .  .2m-l 

(u+2iTj) 


du  . 


A2m-l(0;u) 


From  (10.2)  and  (15.9)  we  have  after  some  easy  calculations 

'n 

(18.7)  < 


2m-l,2(erLU)  =  "  (2m-1)!(2m-2)!(e1U-l)2m  2  H2 


A2m-l(0;elU) 

(2m-l) ! 


/  iu_  .2m 

=  (2m-l)  !  (2m-2)  !  A  .(0;u)  . 

2iu  2m- 1 

e 


It  follows  from  (18.6)  and  (18.7),  after  some  calculations,  that 


(18.8)  N  (x)  = 
o 


(2m-l) ! ( 2m- 2 ) ! (-1) 

2  TT 


m 


iux  ,  .  u  2m 

e  (2  sm  II) 


2m  L  (  0  ..2m- 2 

u  -°°  (u+2ttj) 


00 

—  y 

!m-l  ^ 


2m-l  L  .  ..2m-l 

u  -00  (u+2irj ) 


du, 


from  which  we  obtain  (18.3). 

Similarly  starting  with  (17.16),  we  prove  (18.4).  □ 

Remark.  The  Fourier  transforms  of  the  fundamental  splines  for  C.L.I.P. 
can  also  be  computed  by  the  same  method.  For  instance  if  we  denote  by 
(s)  (s  =  0,2)  the  fundamental  splines  of  (0,2)  interpolation 

2m- 1 , 2 , s 

by  cardinal  g-splines  (r  =  2  ;  §2),  then  using  the  exponential  Euler 


•  9  Sifv 
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g-splines  of  §16  we  obtain  the  following  representations. 


(18.9) 


L2m-l,2,0('x) 


_1_ 

2tt 


xux 


,2  sin  u/2?m  ,  .  , 

( -  >  P2m-A(u)-( 


u 


.  u  /0  2m-  2 
2  sm  /Z 


u 


p2m- 


-2(u)  ] 


[p2m(u)p2m-4(u)  “  [p2m-2(u)]2] 


and 


(18.10) 


L2m-l,2,2(x) 


_1 

2tt 


xux 


0  .  u,  2m  _  .  u.  2  m- 2 

(- - - - )  -(>0  -  ( - )  p2m(u)J 


u 


P2m-2 


du 


u 


(2  sin  U/2)^ 


p2m(u)p2,n-4(u)  '  [p2m-2(u)1 


du 
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Appendix  I 


TABLE  OF  EXAMPLES 


Type  of  C.L.I.P. 

Null  Spaces 

Characteristic  Polynomials 

(0,2)  I.P.  in 

Sg 

5,2 

S8 

o,2 

n8  2(X)  =  7X2+16X+7 
has  2  negative  zeros 

(0,2)  I0P.  in 

^7,2 

Q8 

57,2 

n8  2(X)  =  11X4-104X3-234X2 
-104X+11  has  2  negative 

and  2  positive  zeros 

(0,1,3)  I.P. 

ln  Sf,3 

57,3 

n8  3(X)  =  2(7X2-20X+7) 
has  2  positive  zeros 

(0,3)  I.P. 

“  S7,2 

3y’°(0,l;3,l) 

n  (0,1;3,1;X)  =  3X4-8X3 
/,Z  2 

-153X  -8X+3  has  2  negative 

and  2  positive  zeros 

(0,5)  I.P. 
in  S7>2 

sy;2(o,i;5,i) 

i 

II  2 (0,1 ; 5 , 1 :  )  =  5X4+176X3 
+478X3+176X+5  has  4  nega¬ 
tive  zeros 

(0,2,3)  I.P. 

Sy’°(0,l;2,2) 

ny  3(0,1;2,2:X)  =  3X3+8X+3 

ln  S7,3 

has  2  negative  zeros 

(0,3,4)  I.P. 

57;°3(o’1;3-2)  i 

n  (0,1; 3,2 :X) =  6(X2+5X+1) 

in  S7  o 

has  2  negative  zeros 

(0,1,2)  I.P. 
n3, 2 

S3,2(0,3;k,0)  , 

1 

n3,2(X)  =  50(X2+5X+1) 

in  ^7,3 

- j 

has  2  negative  zeros 

58  - 


' 

Appendix  II 


EXAMPLES  OF  B-SPLINES 


1.  B-Splines  for  (0,3)  -  Interpolation  by  S  „  . 

2 

(1.1)  The  characteristic  polynomial  is  II  (0 ; 3 : A)  =  3A  +14A+3 
eigenvalues  A^  =  (-7+2/l0)/3  ,  ^2  =  (-7-2/l0)/3  . 

(1.2)  The  spline  Nq(x)  =  8(l-x)3  -  2(2-x)3  -  10(l-x)^  +  5(2-x) 

(x^_0)  ,  Nq(x)  =  Nq(-x)  ,  (x<0)  satisfies  Nq(-1)  =  3  , 

N  (0)  =  14  ,  N  (1)  =  3  ,  N(3)(-l)  =  N(3)(0)  =  N(3)(l)  =  0 
o  o  o  o  o 

has  support  in  (-2,2)  . 


(1.3)  The  spline  N^(x)  =  -  y|-  (1-x)3  -  ~  (2-x)3  +  ~  (1-x)^  + 
j2  (2-x)+  ,  (xK>)  ,  N3(x)  =  -N3(-x)  ,  (x<0)  satisfies  N 
N3(0)  =  N3(l)  =  0  ,  N<3)(-1)  =  3  ,  N^3)(0)  =  14  ,  N^3)(l) 
and  has  support  in  (-2,2)  . 


(lo4) 


The  fundamental  function  L  (x) 

o 

L^(v)  =0  V  V  ,  is  given  by 


r 


satisfying  Lq(v)  =  ' 


L  (x) 
o 


/To 

40 


CO 


I  A 

n=-°° 


n 


1 


N  (x-n) 
o 


with 


,  and 


(-1)  = 
=  3  , 


•> 

x 

if  V^0 
if  V=0 


59 
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The  fundamental  function  L^(x) 
Lg(v)  =  0  ,  V  V  ,  is  given  by 


satisfying 


if  V^O 
if  V=0 


L3(x) 


/To 

40 


N^(x-n) 


2.  B-Splines  for  (0,2)  -  Interpolation  by  S® 

J  ,  z 

(2.1)  The  characteristic  polynomial  is  (A)  =  7A2+16A+7  with  eigen- 

3  ,  Z 

values  Y1  =  (-8  +  /l5)/7  ,  Y2  =  (-8  -  /i5)/7  . 

(2.2)  The  spline  M  (x)  =  -3(x+2)2  +  12(x+l)2  -  18x^  +  12(x-l)2  - 

o  +  +  +  + 

3(x-2)2  +  10(x+2)2  +  20 (x+1) 2  -  60x2  +  20(x-l)2  +  10(x-2)2 

+  -r  +  +  +  + 

satisfies  M  (-1)  =  7  ,  M  (0)  =  16  ,  M  (1)  =  7  ,  M^2)(-l)  = 

M^2\o)  =  M^2' (1)  =  0  ,  and  has  support  in  (-2,2)  . 

(2.3)  The  spline  M„  (x)  =  yr  (x+2)2  +  (x+l)“*  -  3x^  +  (x-1)2  +  x-  (x-2)2 

z  z  +  +  +  +  Z  t 

-  ^  (x+2)2  -  1 3 ( x+1 ) 2  -  33x2  -  13(x-l) 2  -  ^  (x-2)2  satisfies 
M^2)(-l)  =  7  ,  M^2)(0)  =  16  ,  M<2)(1)  =  7  ,  M2(-l)  =  M2(0)  =  M2(l) 
=  0  ,  and  has  support  in  (-2,2)  . 


■ 
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(2.4)  The  fundamental  spline  A^(x)  satisfying  Ao(v)  = < 
('2') 

and  y (V)  =0  V  V  ,  is  given  by 

/15 


0  if  V^O 
1  if  V=0 


Vx)  =  ^  J  V  M0(x‘n)  • 

(2) 

The  fundamental  spline  A2(x)  satisfying  7 (V)  = 

and  A 9 (v)  =0  V  v  ,  is  given  by 


0  if  V^O 
[l  if  V=0 


vlS  00  1  n  I 

A2(x)  =  — j  y|  M2(x-n) 


